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^ ; Abstract 

The renormalization group (RG) method is one of the singular perturbation methods 

which is used in search for asymptotic behavior of solutions of differential equations. 
£h '. In this article, time-independent vector fields and time (almost) periodic vector fields 

\Q \ are considered. Theorems on error estimates for approximate solutions, existence of 

approximate invariant manifolds and their stability, inheritance of symmetries from those 
i — i . for the original equation to those for the RG equation, are proved. Further it is proved that 

the RG method unifies traditional singular perturbation methods, such as the averaging 
W ; method, the multiple time scale method, the (hyper-) normal forms theory, the center 

J ■ manifold reduction, the geometric singular perturbation method and the phase reduction. 

A necessary and sufficient condition for the convergence of the infinite order RG equation 

is also investigated. 

> 1 Introduction 

o , 

ON ■ 

ON ■ Differential equations form a fundamental topic in mathematics and its application to nat- 

ural sciences. In particular, perturbation methods occupy an important place in the theory 
of differential equations. Although most of the differential equations can not be solved 
exactly, some of them are close to solvable problems in some sense, so that perturbation 
methods, which provide techniques to handle such class of problems, have been long 
j> ! studied. 

This article deals with a system of ordinary differential equations (ODEs) on a mani- 

5_3 : fold M of the form 

c3 ' dx 

— =sg(t,x,s), xeM, (1.1) 
at 

which is almost periodic in t with appropriate assumptions (see the assumption (A) in 
Sec.2.1), where s e R or C is a small parameter. 

Since s is small, it is natural to try to construct a solution of this system as a power 
series in s of the form 

x = x(t) = x (t) + sxi(t) + s 2 x 2 (t) + ■■■ . (1.2) 
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Substituting Eq.(1.2) into Eq.(l.l) yields a system of ODEs on x , x\,x 2 , - ■ ■ ■ The method 
to construct x(i) in this manner is called the regular perturbation method. 

It is known that if the function g(t,x,s) is analytic in s, the series (1.2) converges 
to an exact solution of (1.1), while if it is not analytic, (1.2) diverges and no longer 
provides an exact solution. However, the problem arising immediately is that one can 
not calculate infinite series like (1.2) in general whether it converges or not, because it 
involves infinitely many ODEs on xq, x\, x%, • ■ • ■ If the series is truncated at a finite-order 
term in e, another problem arises. For example, suppose that Eq.(l.l) admits an exact 
solution x(t) = sin(ef), and that we do not know the exact solution. In this case, the 
regular perturbation method provides a series of the form 

m = st-^(st) 3 + ^(st) 5 + --- . (1.3) 

If truncated, the series becomes a polynomial in t, which diverges as t — » oo although 
the exact solution is periodic in t. Thus, the perturbation method fails to predict qualita- 
tive properties of the exact solution. Methods which handle such a difficulty and provide 
acceptable approximate solutions are called singular perturbation methods. Many singu- 
lar perturbation methods have been proposed so far [1,3,6,38,39,43,44,47,50] and many 
authors reported that some of them produced the same results though procedures to con- 
struct approximate solutions were different from one another [9,38,43,44,47]. 

The renormalization group (RG) method is the relatively new method proposed by 
Chen, Goldenfeld and Oono [8,9], which reduces a problem to a more simple equation 
called the RG equation, based on an idea of the renormalization group in quantum field 
theory. In their papers, it is shown (without a proof) that the RG method unifies conven- 
tional singular perturbation methods such as the multiple time scale method, the boundary 
layer technique, the WKB analysis and so on. 

After their works, many studies of the RG method have been done [10-13,16,21,22,26- 
32,35,36,42,45,46,48,49,53,55-59]. Kunihiro [28,29] interpreted an approximate solu- 
tion obtained by the RG method as the envelope of a family of regular perturbation so- 
lutions. Nozaki et al. [45,46] proposed the proto-RG method to derive the RG equation 
effectively. Ziane [53] , DeVille et al. [32] and Chiba [10] gave error estimates for ap- 
proximate solutions. Chiba [10] defined the higher order RG equation and the RG trans- 
formation to improve error estimates. He also proved that the RG method could provide 
approximate vector fields and approximate invariant manifolds as well as approximate so- 
lutions. Ei, Fujii and Kunihiro [16] applied the RG method to obtain approximate center 
manifolds and their method was rigorously formulated by Chiba [12]. DeVille et al. [32] 
showed that lower order RG equations are equivalent to normal forms of vector fields, 
and this fact was extended to higher order RG equations by Chiba [11]. Applications to 
partial differential equations are appeared in [16,31,42,45,46,55-59]. 

One of the purposes of this paper is to give basic theorems on the RG method extend- 
ing author's previous works [10-12], in which the RG method is discussed for more re- 
stricted problems than Eq.(l.l). At first, definitions of the higher order RG equations for 
Eq.(l.l) are given and properties of them are investigated. It is proved that the RG method 



2 



provides approximate vector fields (Thm.2.5) and approximate solutions (Thm.2.7) along 
with error estimates. Further, it is shown that if the RG equation has a normally hyper- 
bolic invariant manifold N, the original equation (1.1) also has an invariant manifold N E 
which is diffeomorphic to N (Thms.2.9, 2.14). The RG equation proves to have the same 
symmetries (action of Lie groups) as those for the original equation (Thm.2.12). In addi- 
tion, if the original equation is an autonomous system, the RG equation is shown to have 
an additional symmetry (Thm.2.15). These facts imply that the RG equation is easier to 
analyze than the original equation. An illustrative example to verify these theorems is 
also given (Sec. 2. 5). 

The other purpose of this paper is to show that the RG method extends and unifies 
other traditional singular perturbation methods, such as the averaging method (Sec.4.1), 
the multiple time scale method (Sec.4.2), the (hyper-) normal forms theory (Sec.4.3), the 
center manifold reduction (Sec. 3. 2), the geometric singular perturbation method (Sec. 3. 3), 
the phase reduction (Sec. 3. 4), and Kunihiro's method [28,29] based on envelopes (Sec. 4.4) 
A few of these results were partially obtained by many authors [9,38,43,44,47]. The 
present arguments will greatly reveal the relations among these methods. 

Some properties of the infinite order RG equation are also investigated. It is proved 
that the infinite order RG equation converges if and only if the original equation is invari- 
ant under an appropriate torus action (Thm.5.1). This result extends Zung's theorem [54] 
which gives a necessary and sufficient condition for the convergence of normal forms of 
infinite order. The infinite RG equation for a time-dependent linear system proves to be 
convergent (Thm.5.6) and be related to monodromy matrices in Floquet theory. 

Throughout this paper, solutions of differential equations are supposed to be defined 
for all t e R. 



2 Renormalization group method 

In this section, we give the definition of the RG (renormalization group) equation and 
main theorems on the RG method, such as the existence of invariant manifolds and inher- 
itance of symmetries. An illustrative example and comments on symbolic computation 
of the RG equation are also provided. 

2.1 Setting, definitions and basic lemmas 

Let M be an n dimensional manifold and U an open set in M whose closure U is compact. 
Let g(t, % s) be a vector field on U parameterized by t e R and s e C. We consider the 
system of differential equations 

dx 

— = x = sg(t,x,s). (2.1) 
at 

For this system, we make the following assumption. 

(A) The vector field g(t, x, s) is C 1 with respect to time t e R and C°° with respect to 
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x 6 U and s e S, where £ c C is a small neighborhood of the origin. Further, g is an 
almost periodic function with respect to t uniformly in x € U and s e S, the set of whose 
Fourier exponents has no accumulation points on R. 

In general, a function h{t, x) is called almost periodic with respect to t uniformly in 
x € U if the set 

T(h,6) := {t I \\h(t + T,x)-h(t,x)\\ < 6, Vf e R, Vx e t/} c R 

is relatively dense for any 8 > 0; that is, there exists a positive number L such that 
[a, a + L] n r(/z, 5) ^ for all a e R. It is known that an almost periodic function 
is expanded in a Fourier series as hit, x) ~ £ a n (x)e a "' ', (i = V-T), where /l„ e R is 
called a Fourier exponent. See Fink [20] for basic facts on almost periodic functions. 
The condition for Fourier exponents in the above assumption (A) is essentially used to 
prove Lemma 2.1 below. We denote Mod(/z) the smallest additive group of real numbers 
that contains the Fourier exponents X n of an almost periodic function h(t) and call it the 
module of h. 

Let Xfcli ^gkif, x) be the formal Taylor expansion of sg(t, x, e) in e : 

x = £g x (t, x) + e 2 g 2 {t, x) + ■ ■ ■ . (2.2) 

By the assumption (A), we can show that g t (t, x) (i = 1, 2, • • •) are almost periodic func- 
tions with respect to t e R uniformly in x e U such that Mod(g,) c Mod(g). 

Though Eq.(2.1) is mainly considered in this paper, we note here that Eqs.(2.3) and 
(2.5) below are reduced to Eq.(2.1): Consider the system of the form 

x = f{t, x) + sg(t, x, s), (2.3) 

where f(t, •) is a C°° vector field on U and g satisfies the assumption (A). Let <p t be the 
flow of / ; that is, <p t (x ) is a solution of the equation x = f(t, x) whose initial value is x 
at the initial time t = 0. For this system, changing the coordinates by x = (f t (X) provides 

X = g(t,<Pt(X),s) := sg(t,X,s). (2.4) 

We suppose that 

(B) the vector field g satisfies the assumption (A) and there exists an open set W c U 
such that (p t (W) c U and ip t (x) is almost periodic with respect to t uniformly in x 6 W, 
the set of whose Fourier exponents has no accumulation points. 

Under the assumption (B), we can show that the vector field g(t,X,s) in the right hand 
side of Eq.(2.4) satisfies the assumption (A), in which g is replaced by g. Thus Eq.(2.3) 
is reduced to Eq.(2.1) by the transformation jci->I. 
In many applications, Eq.(2.3) is of the form 

x = Fx + sg(x, s) 

= Fx + sg l (x) + s 2 g 2 (x) + ■■■ , xeC", (2.5) 
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where 

(CI) the matrix F is a diagonalizable n xn constant matrix all of whose eigenvalues are 
on the imaginary axis, 

(C2) each g t (x) is a polynomial vector field on C". 

Then, the assumptions (CI) and (C2) imply the assumption (B) because <p t (x) = e Ft x is 
almost periodic. Therefore the coordinate transformation x = e Ft X brings Eq.(2.5) into 
the form of Eq.(2. 1) : X = se' Ft g(e Ft X, s) := sg(t, X, e). In this case, Mod(g) is generated 
by the absolute values of the eigenvalues of F. Note that any equations x = f(x) with C°° 
vector fields / such that /(0) = take the form (2.5) if we put x h-» ex and expand the 
equations in s. 

In what follows, we consider Eq.(2.1) with the assumption (A). We suppose that the 
system (2.1) is defined on an open set U on Euclidean space M = C". However, all 
results to be obtained below can be easily extended to those for a system on an arbitrary 
manifold by taking local coordinates. Let us substitute x = xq + ex\ + s 2 x 2 + • • • into the 
right hand side of Eq.(2.2) and expand it with respect to s. We write the resultant as 

OO OO 

^ e k g k (t, x + ex\ + s 2 x 2 + •••) = X £ k G k (t, x ,x u ■ ■ ■ , x k -\). (2.6) 

k=i k=i 

For instance, Gi,G 2 , G3 and G 4 are given by 

Gi(t,x ) = gi(t,x ), (2.7) 
dg\ 

G 2 (t,x Q ,x l ) = —(t,x )x l +g 2 (t,x ), (2.8) 
ox 

G 3 (t, x , xux 2 ) = *o)A + ~p~(t, *o)*2 + -!r(t, x Q )xi + g 3 (t, x ), (2.9) 

2 ox 1 ox ox 

G^{t,XQ,X\,X 2 ,X-i) = — ^ (t, Xo)Xj + ^ 2 (t,Xo)X\X 2 + -^—(t,xo)xj 

1 d g 2 2 ^82 °§3 

+ 2~dxF^ t,X 1 + ~b~x ( " t,X ^ X2 + ~dx^ t,X °' )Xl + S4(t, x^AO) 

respectively. Note that G ; (z = 1, 2, • • • ) are almost periodic functions with respect to t 
uniformly in x e such that Mod(G,) c Mod(g). With these G,'s, we define the C°° 
maps R h uf : U -> M to be 

R l (y) = \im- f G l (s,y)ds, (2.11) 

r^oo t J 

u ( '\y) = J (.G l (s,y)-R l (y))ds, (2.12) 
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and 



Rfy) = lim - J (GiCj.y, k^GO, ■ • • , wf 1 ^)) - Jg "^T W-*0>))^, (2-13) 
(Gi( 5 ,y, ■ • ■ , ^- 1} (y)) - g W-fcOO - ^(y))^,(2.l4) 

for / = 2, 3, • • • , respectively, where J' denotes the indefinite integral, whose integral 
constants are fixed arbitrarily (see also Remark 2.4 and Section 2.4). 

Lemma 2.1. (i) The maps R t (i = 1, 2, • • • ) are well-defined (i.e. the limits exist), 
(ii) The maps uf'iy) (z = 1, 2, • • • ) are almost periodic functions with respect to t uni- 
formly my £ U such that Mod(w w ) c Mod(g). In particular, u\ are bounded in t e R. 

Proof. We prove the lemma by induction. Since G\(t,y) = gi(t,y) is almost periodic, it 
is expanded in a Fourier series of the form 

gl (t,y)= Yj ^UeR, (2.15) 

A„eMod(gi) 

where A = 0. Clearly R\(y) coincides with a (y). Thus u~'(y) is written as 

u^(y)= f J]a n (y)e a " s ds. (2.16) 

In general, it is known that the primitive function fh(t, y)dt of an uniformly almost pe- 
riodic function h(t, y) is also uniformly almost periodic if the set of Fourier exponents 
of h(t,y) is bounded away from zero (see Fink [20]). Since the set of Fourier exponents 
of gi(t,y) - Ri(y) is bounded away from zero by the assumption (A), u t (y) is almost 
periodic and calculated as 

u\ l) (y) = V — a n {y)e iAnt + (integral constant). (2.17) 

This proves Lemma 2.1 for i = 1. 

Suppose that Lemma 2.1 holds for i = 1,2, • • • ,k - 1. Since Gt(t,x , ■ ■ ■ ,x^_i) and 
u t^(y)> ' " ' > M f* -1) O0 are uniformly almost periodic functions, the composition Gt(t, y, u^iy), • • • , uf~ l \y)) 
is also an uniformly almost periodic function whose module is included in Mod(g) (see 
Fink [20]). Since the sum, the product and the derivative with respect to a parameter 
y of uniformly almost periodic functions are also uniformly almost periodic (see Fink 
[20]), the integrand in Eq.(2.13) is an uniformly almost periodic function, whose module 
is included in Mod(g). The Rt(y) coincides with its Fourier coefficient associated with 
the zero Fourier exponent. By the assumption (A), the set of Fourier exponents of the in- 
tegrand in Eq.(2.13) has no accumulation points. Thus it turns out that the set of Fourier 
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exponents of the integrand in Eq.(2.14) is bounded away from zero. This proves that 
uf^iy) is uniformly almost periodic and the proof of Lemma 2.1 is completed. ■ 

Before introducing the RG equation, we want to explain how it is derived according to 
Chen, Goldenfeld and Oono [8,9]. The reader who is not interested in formal arguments 
can skip the next paragraph and go to Definition 2.2. 

At first, let us try to construct a formal solution of Eq.(2. 1) by the regular perturbation 
method; that is, substitute Eq.(1.2) into Eq.(2.1). Then we obtain a system of ODEs 

( x Q = 0, 
x\ = Gi(t,x ), 



x n — G n (t, xq, ■ ■ ■ , ^n-i), 



Let xo(t) = y e C" be a solution of the zero-th order equation. Then, the first order 
equation is solved as 



Xi(t) 



f 



G l (s,y)ds = R l (y)t + I (G^y) - R x (y)) ds = R x (y)t + u^(y) 



where we decompose x\(t) into the bounded term u t (y) and the divergence term Ri(y)t 
called the secular term. In a similar manner, we solve the equations on x 2 , x 3 , ■ ■ ■ step by 
step. We can show that solutions are expressed as 



n-\ 



x n (f) = uf(y) + 



Zou - ' 
-7j-(y)R«- k (y) 



t + 0(t 2 ), 



(see Chiba [10] for the proof). In this way, we obtain a formal solution of the form 



n-\ 



du ik) 



x(t) ■= x(t,y) =y + J] £ " u t n) W + E ^ \ Rn(y) + Z -QT^ R n-k(y) 



n=\ 



ra=] 



k=l 



t + 0(t 2 ). 



Now we introduce a dummy parameter r e R and replace polynomials V in the above by 
(t - r) 7 . Next, we regard y = y(r) as a function of r to be determined so that we recover 
the formal solution x{t,y): 



x(t,y) = y(r) + £ s n u^{y{r)) + £ s n R n (y(r)) + g JL-(y( T ))R n _ k (y( T ))j (t-r) + 0((t - r) 2 ). 



1 du<-» 

R„(y(T)) + 2_, 

n=\ n=\ \ 

Since x{t,y) has to be independent of the dummy parameter t, we impose the condition 

d 



dr 



x(t,y) = 0, 



which is called the RG condition. This condition provides 

(n) r 00/ n-1 







n=l 



dy 
dt 



71=1 
( 



(») 



\ CO 



dy 



! k=\ 

Thus we see that y(t) has to satisfy the equation dy/dt = ££Lj e^Rkiy), which gives the 
RG equation. Motivated this formal argument, we define the RG equation as follows: 

Definition 2.2. Along with Rj and up, we define the m-th order RG equation for Eq.(2. 1) 
to be 

y = sR x (y) + s 2 R 2 (y) + ■■■ + s m R m (y\ (2.18) 
and the m-th order RG transformation to be 

a ( t m) (y) =y + su^ty) + ■■■ + s^ty). (2.19) 

Domains of Eq.(2.18) and the map or| m) are shown in the next lemma. 

Lemma 2.3. If \e\ is sufficiently small, there exists an open set V = V(s) c U such that 
a^iy) is a diffeomorphism from V into U, and the inverse (o , | m) )" 1 (^) is also an almost 
periodic function with respect to t uniformly in x. 

Proof. Since the vector field g(t, x, s) is C°° with respect to x and s, so is the map a\ m \ 
Since a| m) is close to the identity map if |e| is small, there is an open set V t c U such 
that or, is a diffeomorphism on V t . Since V f 's are e-close to each other and since or| m) 
is almost periodic, the set V := f] t eR Vt is not empty. We can take the subset V c V if 
necessary so that c^ m) (V) c U. 

Next thing to do is to prove that (a^)' 1 is an uniformly almost periodic function. 
Since a^ m) is uniformly almost periodic, the set 



T{a ( r\ 6) = {r I ||aJS(y) - a| m) (y)|| < 6, W e R, Vv e V} 
is relatively dense for any small 5 > 0. For y £ V, put x = or, (y). Then 



(2.20) 



(mjx-l/- (m) / 



, r ||ar[ m) (y)-a^0;)||<L, +T 5 : 

if r e r(or)" ,; ,5), where L, is the Lipschitz constant of the map (a)" 1 ')' 1 ^. Since a 
is almost periodic, we can prove that there exists the number L : = max feR L ; . Now the 
inequality 

(2.22) 



..('»)/ 



(2.21) 



„(») 



„( m h-ii 



(/») 



(a^x)-(c^V«ll<^ 



holds for any small 8 > 0, r e Tiaf 1 ,S) and x e arJ mj (V). This proves that (a) m ') 1 is an 



(m). 



almost periodic function with respect to ? uniformly in x e aj (V) 
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In what follows, we suppose that the m-th order RG equation and the m-th order RG 
transformation are defined on the set V above. Note that the smaller \e\ is, the larger set 
V we may take. 

Remark 2.4. Since the integral constants in Eqs.(2.11) to (2.14) are left undetermined, 
the m-th order RG equations and the m-th order RG transformations are not unique al- 
though Ri(y) is uniquely determined. However, the theorems described below hold for 
any choice of integral constants unless otherwise noted. Good choices of integral con- 
stants simplify the RG equations and it will be studied in Section 2.4. 



2.2 Main theorems 

Now we are in a position to state our main theorems. 

Theorem 2.5. Let or, be the m-th order RG transformation for Eq.(2.1) defined on 
V as Lemma 2.3. If \s\ is sufficiently small, there exists a vector field S(t,y, s) on V 
parameterized by t and s such that 

(i) by changing the coordinates as x = a^iy), Eq.(2.1) is transformed into the system 



y = eRi(y) + s 2 R 2 (y) + ■■■ + s m R m {y) + s m+1 S(t, y, s), 



(2.23) 



(ii) S is an almost periodic function with respect to t uniformly in y e V with Mod(5 ) c 
Mod(g), 

(iii) S (t, y, e) is C 1 with respect to t and C°° with respect to y and e. In particular, S and 
its derivatives are bounded as s — » and t — » oo. 



Proof. The proof is done by simple calculation. By putting x = a " (y), the left hand 
side of Eq.(2.1) is calculated as 

dx d (m) 

li = dt a < 00 



k=i 



k=i 



dt 



m n (k) \ m 

k=l y ) k=\ 



G k (t,y, u) 



(i) 



7=1 



(2.24) 



On the other hand, the right hand side is calculated as 

oo 

sg(t, a ( ; n \y), s) = Yi ^kiU y + su^iy) + s 2 uf\y) + ■■■) 
k=i 

oo 

= YjS k G k (t,y,u ( t l \y),--- ,uf- l) (y)). 
9 



k=\ 



(2.25) 



Thus Eq.(2.1) is transformed into 



y = 



m du (k) X m \ k ~ X du U) 

id + E e *^r ^ E sk My) + E -^o^oo 



i<=i 



/ k=\ 
(k) \~ l 



,;=i 



dy 



k=i 



dy 



7=1 



+ 



(9w 



fc=i 

(A) 



vA 



y V*=i 



fc=] 



«9y 



w + E<- I E ^) 1 1 1 E + Z^wE «"w 



k=i 



dy 



/ &=m+l 



fc=i 
+ 



/c=l 



/ i=m-k+l 



m Q (k) \ 1 00 

id + J] £k -^(y) E ^ G &>y> "r 1 ^). ■ ■ • . «f _1) 00). 



(2.26) 



/ k=m+\ 



The last two terms above are of order 0(s m+i ) and almost periodic functions because they 
consist of almost periodic functions uf and G ; . This proves Theorem 2.5. ■ 

Remark 2.6. To prove Theorem 2.5 (i),(iii), we do not need the assumption of almost 
periodicity for git, x, s) as long as Rfy) are well-defined and g, w, !) and their derivatives 
are bounded in t so that the last two terms in Eq.(2.26) are bounded. In Chiba [10], The- 
orem 2.5 (i) and (iii) for m = 1 are proved without the assumption (A) but assumptions 
on boundedness of g, u[ and their derivatives. 

Thm.2.5 (iii) implies that we can use the m-th order RG equation to construct ap- 
proximate solutions of Eq.(2.1). Indeed, a curve a t (y(t)), a solution of the RG equation 
transformed by the RG transformation, gives an approximate solution of Eq.(2.1). 

Theorem 2.7 (Error estimate). Let y(t) be a solution of the m-th order RG equation 
and or, the m-th order RG transformation. There exist positive constants sq, C and T 
such that a solution x(t) of Eq.(2.1) with jc(0) = a ™ (y(Q)) satisfies the inequality 

\\x(t)-a ( r\y(t))\\<C\sr, (2.27) 
as long as \s\ < s Q , y(t) e V and < t < T l\e\. 

Remark 2.8. Since the velocity of yit) is of order 0(s), y(0) e V implies y(t) e V for 
< t < T l\e\ unless y(0) is e-close to the boundary of V. If we define so that the 
indefinite integrals in Eqs.(2.12, 14) are replaced by the definite integrals J , a" is the 
identity and a ( ™\y(0)) = y(0). 

Proof of Thm.2.7. Since a| m) is a diffeomorphism on V and bounded in t e R, it is 
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sufficient to prove that a solution y(t) of Eq.(2.18) and a solution y{t) of Eq.(2.23) with 
y(0) = y(0) satisfy the inequality 



HKO-y(OII<C|er, 0<t<T/\e\, (2.28) 

for some positive constant C. 

Let Li > be the Lipschitz constant of the function Ri(y) + sRiiy) + • • • + e!"~ x R m (y) 
on V and L 2 > a constant such that sup, eR ve y \\S(t,y, e)\\ < L 2 . Then, by Eq.(2.18) and 
Eq.(2.23), y(t) and y(t) prove to satisfy 

HK0-y(0ll<^i f \m-y(s)\\ds + L 2 e m+i t. (2.29) 
Jo 

Now the Gronwall inequality proves that 

\\m-y(t)\\<^"\e ELlt -l). (2.30) 

The right hand side is of order 0(e m ) ifO<t<T/e. ■ 

In the same way as this proof, we can show that if Ri(y) = ■ ■ ■ = Rk(y) = holds with 
k < m, the inequality (2.27) holds for the longer time interval < t < T l\e\ k+l . This fact 
is proved by Murdock and Wang [41] for the case k = 1 in terms of the multiple time 
scale method. 

We can also detect existence of invariant manifolds. Note that introducing the new 
variable s, we can rewrite Eq.(2.1) as the autonomous system 

dx 

— = sg(s, x, e), 
dt 

Then we say that Eq.(2.31) is defined on the (s, x) space. 

Theorem 2.9 (Existence of invariant manifolds). Suppose that R\(y) = ■■■ = Rk-i(y) = 

and s k Rt(y) is the first non-zero term in the RG equation for Eq.(2.1). If the vector field 
Rk(y) has a boundaryless compact normally hyperbolic invariant manifold N, then for 
sufficiently small e > 0, Eq.(2.3 1) has an invariant manifold A^ e on the (s, x) space which 
is diffeomorphic to R x N. In particular, the stability of Af e coincides with that of N. 

To prove this theorem, we need Fenichel's theorem : 
Theorem (Fenichel [18]). Let M be a C 1 manifold and X(M) the set of C 1 vector fields 
on M with the C 1 topology. Suppose that / e X(M) has a boundaryless compact normally 
hyperbolic /-invariant manifold N c M. Then, the following holds: 

(i) There is a neighborhood tl c X(M) of / such that there exists a normally hyperbolic 
g-invariant manifold N g c M for any g eU. The N g is diffeomorphic to N. 
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(ii) If ||/ - g| | ~ 0(e), N g lies within an 0(e) neighborhood of N uniquely. 

(iii) The stability of N E coincides with that of N. 

Note that for the case of a compact normally hyperbolic invariant manifold with 
boundary, Fenichel's theorem is modified as follows : If a vector field / has a compact 
normally hyperbolic invariant manifold N with boundary, then a vector field g, which is 
C 1 close to /, has a locally invariant manifold N g which is diffeomorphic to N. In this 
case, an orbit of the flow of g on N g may go out from N g through its boundary. Accord- 
ing to this theorem, Thm.2.9 has to be modified so that _/V e is locally invariant if N has 
boundary. 

See [18,24,51] for the proof of Fenichel's theorem and the definition of normal hy- 
perbolicity. 

Proof of Thm.2.9. Changing the time scale as 1 1-» t/s k and introducing the new variable 
s, we rewrite the k-th order RG equation as 



ds 

di = ' 



(2.32) 



and Eq.(2.23) as 

dy 



— = R k (y) + sR k+i (y) + 
dt 
ds 

di = ' 



+ s m - k R m (y) + £ m+i ~ k S(s/s k ,y,s), 



m+l-k i 



(2.33) 



respectively. Suppose that m > 2k. Since S is bounded in s and since 

^-s m+l - k S( S /s k ,y,s) ~ 0(s k+l ), ^-s m+l - k S(s/s k ,y,s) - 0(e), 
oy as 



(2.34) 



Eq.(2.33) is e-close to Eq.(2.32) on the (s,y) space in the C 1 topology. 

By the assumption, Eq.(2.32) has a normally hyperbolic invariant manifold RxA^on 
the (s,y) space. At this time, Fenichel's theorem is not applicable because R x Af is not 
compact. To handle this difficulty, we do as follows: 

Since S is almost periodic, the set 



T(S,6) := {t\\\S((s - r)le k ,y,s)- S(s/e k ,y,s)\\ < 5, VseR} 



(2.35) 



is relatively dense for any small 8 > 0. Let us fix 5 so that it is sufficiently smaller than 
e and fix r e T(S,6) arbitrarily. Then W := [0, r] x N is a compact locally invariant 
manifold of Eq.(2.32) with boundaries {0} x Af and {t} x Af (see Fig.l). 

Now Fenichel's theorem proves that Eq.(2.33) has a locally invariant manifold W £ 
which is diffeomorphic to W and lies within an 0(e) neighborhood of W uniquely. 
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To extend W s along the s axis, consider the system 



( 



y = R k (y) + sR k+1 (y) + ■■■ + s m ' k R m {y) + s' 



S((s-r)/s k ,y,s), 



(2.36) 



s = 1. 



Since the above system is 5-close to Eq.(2.33), it has a locally invariant manifold W e j, 
which is diffeomorphic to W E . By putting s = s - r, Eq.(2.36) is rewritten as 



is a locally invariant manifold of Eq.(2.33). Since W St g is 5-close to W E and since 6 <sc s, 
both of W £ n {s = t} and K n {s = r} are e-close to W. Since an invariant manifold of 
Eq.(2.33) which lies within an 0(e) neighborhood of W is unique by Fenichel's theorem, 
K P\{s = t} has to coincide with W E C\{s = t}. This proves that K is connected to W s and 
KU W E gives a locally invariant manifold of Eq.(2.33). 

This procedure is done for any t g T(S,6). Thus it turns out that W E is extended 
along the s axis and it gives an invariant manifold N E - R x N of Eq.(2.33). An invariant 
manifold A^ £ of Eq.(2.1) is obtained by transforming A^ £ by a ( " r \ 

Note that by the construction, projections of the sets N £ C\ {s = r}, r 6 T(S,6) on to 
the y space are 5-close to each other. This fact is used to prove the next corollary. ■ 

The next corollary (ii) and (iii) for k = 1 are proved in Bogoliubov, Mitropolsky [6] 
and Fink [20] and immediately follow from Thm.2.9. 

Corollary 2.10. Suppose that R\(y) = ■ • ■ = Rk-i(y) = and e*i?^(y) is the first non-zero 
term in the RG equation for Eq.(2.1). For sufficiently small e > 0, 

(i) if the vector field Rtiy) has a hyperbolic periodic orbit ja(t), then Eq.(2.1) has an 
almost periodic solution with the same stability as jo(t), 

(ii) if the vector field R k (y) has a hyperbolic fixed point y , then Eq.(2.1) has an almost 
periodic solution y £ (i) with the same stability as jq such that Mod(y £ ) c Mod(g), 

(iii) if the vector field Ruiy) has a hyperbolic fixed point y and if g is periodic in t with a 
period T, then Eq.(2.1) has a periodic solution y £ (f) with the same stability as y and the 
period T (it need not be the least period). 

Proof. If Rkiy) has a periodic orbit, Eq.(2.33) has an invariant cylinder A^ £ on the (s,y) 
space as is represented in Fig.l. To prove Corollary 2.10 (i), at first we suppose that 
g(t, x, s) is periodic with a period T. In this case, since S(t,y, s) is a periodic function 
with the period T, N £ is periodic along the s axis in the sense that the projections S 1 : = 
N E P\{s = mT] give the same circle for all integers m. Let y = y(t), s = t be a solution of 
Eq.(2.33) on the cylinder. Then y(mT), m = 0, 1, • • • gives a discrete dynamics on S 1 . If 
y(mT) converges to a fixed point or a periodic orbit as m — » oo, then y(t) converges to a 




.m+l-k 



S(s/e k ,y,s), 



(2.37) 



and it takes the same form as Eq.(2.33). This means that the set 



K:={( S ,y)\( S -r,y)eW £ J 
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Fig. 1 : A schematic view of the proof for the case that N is a circle. The W E is e-close to 
W and W StS is enclose to W e . The K is the "copy" of W s>6 . 



periodic function as t —* oo. Otherwise, the orbit of y(mT) is dense on S 1 and in this case 
y(t) is an almost periodic function. A solution of Eq.(2.1) is obtained by transforming 
y(t) by the almost periodic map a* m) . This proves (i) of Corollary 2.10 for the case that g 
is periodic. 

If g is almost periodic, the sets N E C\{s = t} give circles for any r 6 T(S,6) and they are 
<5-close to each other as is mentioned in the end of the proof of Thm.2.9. In this case, there 
exists a coordinate transformation Y = (p(y, t) such that the cylinder A^ £ is straightened 
along the s axis. The function tp is almost periodic in t because \\ipiy, t + r) - (p(y, Oil is of 
order 0(6) for any r e T(S,6). Now the proof is reduced to the case that g is periodic. 

The proofs of (ii) and (iii) of Corollary 2.10 are done in the same way as (i), details 
of which are left to the reader. ■ 

Remark 2.11. Suppose that the first order RG equation sRi(y) + does not have nor- 
mally hyperbolic invariant manifolds but the second order RG equation sR\(y) + £ 2 R 2 (y) 
does. Then can we conclude that the original system (2.1) has an invariant manifold with 
the same stability as that of the second order RG equation? Unfortunately, it is not true in 
general. For example, suppose that the RG equation for some system is a linear equation 
of the form 



y/s 



1\ I 0\ 2 0. 

OP"* l y + £ 4 0^ + -'-' ^ eR "- 



(2.38) 



The origin is a fixed point of this system, however, the first term has zero eigenvalues and 
we can not determine the stability up to the first order RG equation. If we calculate up to 
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the second order, the eigenvalues of the matrix 



° ^ 41 ? (2.39) 



o o; \o i 

are —s (double root), so that y = is a stable fixed point of the second order RG equation 
if s > 0. Unlike Corollary 2.10 (ii), this does not prove that the original system has a 
stable almost periodic solution. Indeed, if we calculate the third order RG equation, the 
eigenvalues of the matrix in the right hand side of Eq.(2.38) are 3s and —s. Therefore the 
origin is an unstable fixed point of the third order RG equation. This example shows that 
if we truncate higher order terms of the RG equation, stability of an invariant manifold 
may change and we can not use sR\(y) + s 2 R 2 (y) to investigate stability of an invariant 
manifold as long as R\(y) ± 0. This is because Fenichel's theorem does not hold if the 
vector field / in his theorem depends on the parameter s. 

Theorems 2.7 and 2.9 mean that the RG equation is useful to understand the properties 
of the flow of the system (2.1). Since the RG equation is an autonomous system while 
Eq.(2. 1) is not, it seems that the RG equation is easier to analyze than the original system 
(2.1). Actually, we can show that the RG equation does not lose symmetries the system 
(2.1) has. 

Recall that integral constants in Eqs.(2.12, 14) are left undetermined and they can 
depend on y (see Remark 2.4). To express the integral constants Bfy) in Eqs.(2.12, 14) 
explicitly, we rewrite them as 

u ( t 1 \y) = B l (y) + j (G l ( S ,y)-R l (y))ds, 



and 



(g,(5, y, u?(y), • • • , « ( ;'- 1} Cy)) - -^-(y)Ri-k(y) - Ri(y))ds, 

for i = 2, 3, • • • , where integral constants of the indefinite integrals in the above formulas 
are chosen to be zero. 

Theorem 2.12 (Inheritance of symmetries). Suppose that an e-independent Lie group 
H acts on U c M. If the vector field g and integral constants B t (y), i = 1, • • • ,m - 1 in 
Eqs.(2.12, 14) are invariant under the action of H; that is, they satisfy 

g(t, hy, s) = ^(y)g(t, y, s), BtQiy) = J(y)^(y), (2.40) 
ay oy 

for any h e H,y e U, t e R and e, then the m-th order RG equation for Eq.(2.1) is also 
invariant under the action of H. 

Proof. Since h e H is independent of s, Eq.(2.40) implies 

gi(t,hy) = ^(y) gi (t,y), (2.41) 
dy 
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for / = 1, 2, • • • . We prove by induction that Rfy) and Uj(y), i = 1, 2, ■ ■ ■ , are invariant 
under the action of H. At first, R\{hy), h e H is calculated as 

Riihy) = lim- f G l {t,hy)ds 
• t J 

I C'dh _ 5/i 

dy 

Next, is calculated in a similar way: 



= lim 7 f ^(y) 
f^°° t J Oy 



)Bi(y) + —(y) J (GjU^-i?^))^ 



5/i 
<9y 



dy 

Suppose that 7?jt and are invariant under the action of H for k = 1, 2, • • • , i - 1. Then, 
it is easy to verify that 

dwf , dh 3uf (dh ^ 



(hy) = —(y)^-(y)\—(y)\ , (2.43) 
oy oy Oy \0y ) 

dh 

G k (hy, uf\hy), • • • , uf~\hy)) = ^-(y)G k (y, u { ?(y), ■■■ , uf~ l) (y)), (2.44) 

for k = 1,2, 1. These equalities and Eqs.(2.13), (2.14) prove Theorem 2.12 by a 

similar calculation to Eq.(2.42). ■ 

2.3 Main theorems for autonomous systems 

In this subsection, we consider an autonomous system of the form 

x = f(x) + sg(x,s) 

= f{x) + egi(x) + s 2 g 2 {x) + ■■■ , x 6 U c M, (2.45) 

where the flow ip, of / is assumed to be almost periodic due to the assumption (B) so that 
Eq.(2.45) is transformed into the system of the form of (2.1). For this system, we restate 
definitions and theorems obtained so far in the present notation for convenience. We also 
show a few additional theorems. 

Definition 2.13. Let <p t be the flow of the vector field /. For Eq.(2.45), define the C°° 
maps R h hf : U -> M to be 

Ri(y) = ]xm\ f (DtpX 1 Gi(s,<p s (y))ds, (2.46) 
hl l) (y) = (D<p t ) y j ((Dip^Giis, (p s (y)) - R&j) ds, (2.47) 
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and 



Rfy) = lim - [((Dcp^Gtis, <p s (y), h™(y), • • • , h^(y)) 

i-\ 

-{Dcp s ); l Y J i.Dhf) y R i . k {y))ds, (2.48) 

k=\ 

hf(y) = (D<p t ) y f t ((D<p s r y 1 G i (s,<p s (y),h?(y),--- ,h^(y)) 

i-i 

-(Dtps); 1 Y^DhfXR^iy) - Rfy))ds, (2.49) 

for i = 2, 3, • • • , respectively, where (Dhf^) y is the derivative of hf\y) with respect to 
y, (D(f t ) y is the derivative of (p t (y) with respect to y, and where G, are defined through 
Eq.(2.6). With these R, and h t , define the m-th order RG equation for Eq.(2.45) to be 

y = eRiiy) + e 2 R 2 (y) + ■■■ + e m R m (y), (2.50) 

and define the m-th order RG transformation to be 

a ( t n) (y) = <P,(y) + ^(y) + ■■■ + s m h^{y), (2.51) 

respectively. 

In the present notation, Theorems 2.5 and 2.7 are true though the relation Mod(S) c 
Mod(g) in Thm.2.5 (ii) is replaced by Mod(S) c Mod(^,). Note that even if Eq.(2.45) is 
autonomous, the function S depends on t as long as the flow ip t depends on t. 

Theorem 2.9 is refined as follows: 

Theorem 2.14 (Existence of invariant manifolds). Suppose that R\(y) = ••• = 
Rk-i(y) = an d s k Rk(y) is the first non-zero term in the RG equation for Eq.(2.45). 
If the vector field R^iy) has a boundaryless compact normally hyperbolic invariant man- 
ifold N, then for sufficiently small s > 0, Eq.(2.45) has an invariant manifold iV £ , which 
is diffeomorphic to N. In particular, the stability of Af e coincides with that of N. 

Note that unlike Thm.2.9, we need not prepare the (s, x) space, and the invariant 
manifold Af e lies on M not R x M. This theorem immediately follows from the proof of 
Thm.2.9. Indeed, Eq.(2.33) has an invariant manifold A^ £ - RxA^on the (s,y) space 
as is shown in the proof of Thm.2.9. An invariant manifold of Eq.(2.45) on the (s, x) is 
obtained by transforming N £ by the RG transformation. However, it has to be straight 
along the s axis because Eq.(2.45) is autonomous. Thus its projection onto the x space 
gives the invariant manifold iV e of Eq.(2.45) (see Fig. 2). 

For the case of Eq.(2.1), the RG equation is simpler than the original system (2.1) 
in the sense that it has the same symmetries as (2.1) and further it is an autonomous 
system while (2.1) is not. In the present situation of (2.45), Theorem 2.12 of inheritance 
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flow of the RG equation 



flowofEq.(2.33) 



flow of Eq.(2.45) 



Fig. 2: A schematic view of the proof for the case that N is a circle. The projection of 
the straight cylinder on to the x space gives an invariant manifold of Eq.(2.45). 

of symmetries still holds as long as the assumption for g is replaced as "the vector field 
/ and g are invariant under the action of a Lie group H". However, since Eq.(2.45) is 
originally an autonomous system, it is not clear that the RG equation for Eq.(2.45) is 
easier to analyze than the original system (2.45). The next theorem shows that the RG 
equation for Eq.(2.45) has larger symmetries than Eq.(2.45). 

To express the integral constants Bfy) in Eqs.(2.47, 49) explicitly, we rewrite them as 

h ( ?(y) = (P<p t ) y B x (y) + (Dcp t ) y j ((D^G^s, <p s (y)) - R^y)) ds, 

and 

hf(y) = (Pip t ) y Bfy) + (Dcp t ) y JliD^X'Gtis, <p s (y), h?(y), h^(y)) 

-{Dtp s ); x Yp h f)y R i-M - Ri(y))ds, 
k=i 

for i = 2,3,'", where integral constants of the indefinite integrals in the above formulas 
are chosen to be zero. 

Theorem 2.15 (Additional symmetry). Let (p t be the flow of the vector field / defined 
on U c M. If the integral constants B t in Eqs.(2.47, 49) are chosen so that they are 
invariant under the action of the one-parameter group {tp, : U —* M\ t e R}, then the RG 
equation for Eq.(2.45) is also invariant under the action of the group. In other words, R { 
satisfies the equality 

Ri(<p t (y)) = (D(p t )yRi(y), (2.52) 

for? = 1,2,-". 

Proof. Since Eq.(2.45) is autonomous, the function Gu{t, xq,-- , x k -\) defined through 
Eq.(2.6) is independent of t and we write it as Gk(x , ■■■ , x^-i). We prove by induction 
that equalities Ri(<p t (y)) = (Dip^Rfy) and hf(ip f (y)) = hf +f (y) hold for i = 1, 2, • • • . For 
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all s' 6 R, R\(ip S '(y)) takes the form 

1 r' 

Ri((p s >(y)) = lim- (Dtp s ) l() G\((p s o (p s ,(y))ds 

= (D<p s ,) y \im- \ (D<p s+s ,)~ 1 G 1 ((p s+s ,(y))ds. 

■ t^oo t J 

Putting s + s' = s", we verify that 

1 r t+s ' 

RMAy)) = (Dip s ,) y ]im- {Dip^G^iyW 

■ t^oo t J 

i r t+s ' 

= (DtpMW + iDtp^Mm- {Dip^fG^AyW 
= (DipMiy). 

The hr t ((p S '(y)) is calculated in a similar way as 

= (£W) y fliO0 + (D<p t+S ,) y J ((D^X'G^^iy)) - Rt(y)) ds. 
Putting s + s' = s" provides 

/t+s' 
((Dtp^GMAyV-Riiy^ds" 

= h%(y). (2.53) 

Suppose that R k (<p t (y)) = (D<p t ) y R k (y) and hf(cp f (y)) = (y) hold for k = 1, 2, ■ ■ ■ , i - 1 . 
Then, Rj((p s >(y)) is calculated as 

Rt(<p,(y)) = lim - f ((Dip s o ^,(y), ^(yvOO), ■ • • , fc^fovGO)) 

i-l 

= (D^),limy J((D^); 1 G I -(^(y),^ ) y (y),.-. 

(-1 

-(D(p s+S , )~ 1 ^ (Dhfl s , ) y Ri-k(y))d s . 



k=i 
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Putting s + s' = s" provides 

Rt^Ay)) = (D<p s ,) y liml f ((D<p s ,,);%(<p s ,,(y),h%) r -- ,h ( <; l) (y)) 

1-1 

-(D(p s » )~ 1 ^ (Dh^)yRi-ic(y))d s" 
k=i 

= (D<p s ,)yRi(y) + (Dip,),, lim - f ((Dtp^G^iy), h%), ■■■ , h^(y)) 

1-1 

= (Dif^yRfy). 

We can show the equality hy((p t >(y)) = h® f (y) in a similar way. ■ 

2.4 Simplified RG equation 

Recall that the definitions of i?, and uf given in Eqs.(2.1 1) to (2.14) include the indefinite 
integrals and we have left the integral constants undetermined (see Rem. 2. 4). In this sub- 
section, we investigate how different choices of the integral constants change the forms 
of RG equations and RG transformations. 

Let us fix definitions of R, and uy, i = 1,2,--- by fixing integral constants in Eqs. (2. 12, 
14) arbitrarily (note that Rj and up are independent of the integral constants in Eqs. (2.1 1, 
13)). For these Rj and uf , we define Rj and uy to be 

^00 = lim- f G l (s,y)ds, (2.54) 

rf t l) (y) = B 1 (y) + J (G l (s,y)-R l (y))ds, (2.55) 

and 

Rfy) = lim - (Gi(s,y,ltP(y), ■ ■ ■ , «f 1} (y)) - J] -^(y)Ri-k(y))ds, (2.56) 

r^oo t J ^ ay 

(Gfay^w, ■ ■ ■ rft 1} (y)) - 2j -^(y)Ri-k(y)-Ri(y))W) 

k=i y 

for i = 2, 3, • • •, respectively, where integral constants in the definitions of w; are the 
same as those of up and where B h i = 1, 2, • • • are arbitrary vector fields which imply 
other choice of integral constants. Along with these functions, We define the m-th order 
RG equation and the m-th order RG transformation to be 

y = sR i (y) + --- + s m R m (y), (2.58) 



20 



and 

5j»)(y) = y + £ ^\ y) + ... + s "rf n \y), (2-59) 

respectively. Now we have two pairs of RG equations-transformations, Eqs.(2.18),(2.19) 
and Eqs. (l2.58l) . (l2.59l) . Main theorems described so far hold for both of them for any 
choices of 5,(y)'s except to Thms. 2.11 and 2.15, in which we need additional assump- 
tions for 5,(y)'s as was stated. 

Let us examine relations between R it and Rj, u{ . Clearly Ri coincides with R x . 
Thus u> t ' is given as ir t '(y) = u, (y) + Bi(y). According to the definition of G 2 (Eq.(2.8)), 
i?2 is calculated as 

R 2 (y) = lim I f y)^(y) + g 2 (j, y) - -jMiGo) 

r^oo t J \ oy oy ) 

1 C l i d d ^ \ 

= lim 7 -^(j.yV.^ + ^.y)--??-^)^ 

f ^°° ? J \ oy oy I 

+ lim- ^.(^y^iCy)-— i(y)/?i(y) Wis 
f^co ? J \ ay ay / 

= R 2 (y) + 1 r 1 (y)B l (y)-—±(y)Ri(y)- (2.60) 
oy oy 

If we define the commutator [ • , • ] of vector fields to be 

dBi dRi 
[B u Ri\(y) = -r 1 (y)R l (y) - -^ 1 (y)B l (y), (2.61) 

oy oy 

Eq. (12.60l) is rewritten as 

Ri(y) = R2(y)-[BuRi](y). (2.62) 

Similar calculation proves that R h i = 2, 3, • • • are expressed as 

Ri(y) = Ri(y) + ,Ri-i,B u -- ,5,_ 2 )(y) - [^/^(y), (2.63) 

where P,- is a function of • • • i and B\, ■ ■ ■ ,5,-2- See Chiba[ll] for the proof. 
Thus appropriate choices of vector fields 5,(y), i = 1, 2, • • • may simplify Rj, i = 2, 3, • • • 
through Eq. (l2331 . 

Suppose that i = 1,2, •• • are elements of some finite dimensional vector space 
V and the commutator [ • ,R X ] defines the linear map on V. For example if the function 
g(t, x, s) in Eq.(2.1) is polynomial in x, V is the space of polynomial vector fields. If 
Eq.(2.1) is an ^-dimensional linear equation, then V is the space of all n x n constant 
matrices. Let us take a complementary subspace C to Im [ • ,R\] into V arbitrarily: V = 
Im[-,7?i] 0C. Then, there exist B { e V, i = 1,2, • • • ,m - 1 such that R { e C for 
i = 2, 3, • • • , m because of Eq. (l2.63l) . If Rj e C for / = 2, 3, • • • , m, we call Eq. (l2.58l) the 
m-th order simplified RG equation. See Chiba[l 1] for explicit forms of the simplified RG 
equation for the cases that Eq.(2.1) is polynomial in x or a linear system. In particular, 



21 



they are quite related to the simplified normal forms theory (hyper-normal forms theory) 
[2,39,40]. See also Section 4.3. 

Since integral constants 5,'s in Eqs. (l2.55l) and (12.571) are independent of t, we can 
show the next claim, which will be used to prove Thm.5.1. 

Claim 2.16. RG equations and RG transformations are not unique in general be- 
cause of undetermined integral constants in Eqs.(2.12) and (2.14). Let a\ m) and c^ m) 
be two different RG transformations for a given system (2.1). Then, there exists a 
time-independent transformation <p(y, s), which is C°° with respect to y and e, such that 
00 = a t ° 0Cy, s). Conversely, af^ o <p(y, s) gives one of the RG transformations for 
any C°° maps <p(y, s). 

Note that the map (p(y, s) is independent of t because it brings one of the RG equations 
into the other RG equation, both of which are autonomous systems. According to Claim 
2.16, one of the simplest way to achieve simplified RG equations is as follows: At first, 
we calculate 7?, and u\ l) by fixing integral constants arbitrarily and obtain the RG equation. 
It may be convenient in practice to choose zeros as integral constants (see Prop.2.18 
below). Then, any other RG equations are given by transforming the present RG equation 
by time-independent C°° maps. 

2.5 An example 

In this subsection, we give an example to verify the main theorems. See Chiba[10,l 1] for 
more examples. Consider the system on R 2 



where e > 0, k > and a> > are parameters. Changing the coordinates by (X l ,X 2 ) = 
(sxi,sx 2 ) yields 



X2 = —X\ + £(*2 _ ^1^2) + e X%, 

Diagonalizing the unperturbed term by introducing the complex variable z as x\ = z + 
Z, x 2 = i(z - z) may simplify our calculation : 



where i = V— L Let us calculate the RG equation for the system (12.651) or (12.661) . In this 
example, all integral constants in Eqs.(2.47, 49) are chosen to be zero. 



( 



X\ = X 2 +Xj + s 2 k sm(u)t), 

X2 — — X\ H" S X2 — X1X2 "I - 1 



(2.64) 




z = iz+- (i(z - zf - 2z 2 + 2zz + k sin(wf)) + — (z - z), 

£ S 2 

z = -iz + - {-i(z - zf - 2z 2 + 2zz + k sin(wf)) - — (z - z), 



(2.66) 
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(i) When co ± 1, 2, the second order RG equation for Eq. (l2.66l) is given as 



1 



h = 2 sl{yx ~ 3y2iy2 



16/ 



3 -yiyi), 

h = 2 s 0^2 - 3yiy 2 + — yiy 2 ), 



(2.67) 



where the first order RG equation Ri vanishes. Note that it is independent of the time 
periodic external force k s'm(cot). Thus this RG equation coincides with that of the au- 
tonomous system obtained by putting k = in Eq. (l2.66l) and Theorem 2.15 is applicable 
to Eq. (l2.67l) . Indeed, since the above RG equation is invariant under the rotation group 



Oi,)>2) ^ (e' T y\,e lT y 2 ), putting y x 



if) 

re', y 2 



1 



re results in 



r = -s 2 r(l - 3r 2 ), 



(2.68) 



and it is easily solved. We can verify that this RG equation has a stable periodic orbit 
r = VI /3 if s > 0. Now Corollary 2.10 (i) proves that the original system (12.651) has a 
stable almost periodic solution if £ > is small (see Fig. 3). If k = and Eq. (l2.65l) is 
autonomous, then Thm.2.14 is applied to conclude that Eq. (l2.65l) has a stable periodic 
orbit. 



Xi o 




Fig. 3: Numerical results of the system (12.651) and its RG equation (12.681) for co = 3, k = 
1.8 and s = 0.01. The red curve denotes the stable periodic orbit of the RG equation and 
the black curve denotes the almost periodic solution of Eq. d2.65l) . They almost overlap 
with one another. 
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(ii) When to = 2, we prefer Eq. (l2.65l) to Eq. (l2.66l) to avoid complex numbers. The second 
order RG equation for Eq. (l2.65l) is given by 



\ yi = — illy, - 9y\ - I6y]y 2 - 9y x y\ -I6y 3 2 - k(6 yi + 4y 2 j) , 

I 2 f (2.69) 

[ y 2 = ^ (l2y 2 - 9y\ + I6 yi y 2 2 - 9y\y 2 + \6y\ - k(4 yi - 6y 2 )) . 

Since it depends on the time periodic term ks'm(a) 2 t), Thm.2.15 is no longer applicable 
and to analyze this RG equation is rather difficult. However, numerical simulation shows 
that Eq. fl2.69l) undergoes a typical homoclinic bifurcations (see Chow, Li and Wang [14]) 
whose phase portraits are given as Fig.4. 




£-1.8 



Fig. 4: Phase portraits of the RG equation (|2.69b . 

Let us consider the case k = 1.8. In this case, the RG equation has one stable periodic 
orbit and two stable fixed points. Thus Corollary 2.10 proves that the original system 
(12.651) has one almost periodic solution y(t) and two periodic solutions with the period n 
(see Fig. 5; in the arXiv version, Fig. 5 is omitted because of the limitation of size). 

(iii) When a> = 1, the second RG equation for Eq. (l2.65l) is given by 

yi = — (l2 yi - 9y] - \6y\y 2 - 9y,y\ - 16$ , 

Ik s 2 (2 - 70) 

J 2 = T + h ( 12y2 _ 9y2 + l6yiy2 ~ 9y2yy2 + l6y ^ ■ 

In this case, the first order RG equation does not vanish. For small s, Eq. fl2.70l) has a 
stable fixed point y = y (k, s), which tends to infinity as e — » 0. For example if k = 1.8 
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ande = 0.01, y is given by y ~ (-4.35, 2.31) (see Fig. 6). Therefore, the original system 
has a stable periodic orbit whose radius tends to infinity as e —* 0. 

5 
4 

3 
2 
1 

x 2 
-1 
-2 
-3 
-4 
■5 

■5 -4 -3 -2 -1 1 2 3 4 5 6 

Fig. 6: Numerical results of the system (12.651) and its RG equation (12.701) for a> = 1, k = 
1.8 and s = 0.01. The red cross point denotes the stable fixed point of the RG equation. 
The black curve denotes the periodic solution of Eq. (l2.65l) . 



2.6 A few remarks on symbolic computation 

One can calculate the RG equation and the RG transformation by using symbolic com- 
putation softwares, such as Mathematica and Maple, with formulas (2.11) to (2.14). In 
this subsection, we provide a few remarks which may be convenient for symbolic com- 
putation. 

It is not comfortable to compute the limits in formulas (2.11) and (2.13) directly by 
symbolic computation softwares because it takes too much time. The next proposition is 
useful to compute them. 

Proposition 2.17. Suppose that F(t, y) and its primitive function are almost periodic 
functions in t. Then, lim,^ j J F(t, y)dt gives a coefficient of a linear term of J F(t, y)dt 
with respect to t. 

This proposition is easily proved because F is expanded in a Fourier series as F(t, y) ~ 
X a n (y)e Unt . Thus, to obtain 2?,(y)'s, we compute integrals in Eqs. (2.1 1,13) and extract 
linear terms with respect to t. To do so, for example in Mathematica, the command 
Coefficient [Integrate [F [t , y] , t] , t] is available, where F [t , y] is the integrand 
in Eqs.(2.11, 13). 
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When computing the integrals in Eqs.(2.11) to (2.14), we recommend that the inte- 
grands are expressed by exponential functions with respect to t such as e lAt if they include 
trigonometric functions such as sin At, cos At (in Mathematica, it is done by using the 
command TrigToExp). It is because if they include trigonometric functions, softwares 
may choose unexpected integral constants while if they consist of exponential functions, 
then zeros are chosen as integral constants. If all integral constants in Eqs.(2.12), (2.14) 
are zeros, the second term in the integrand in Eq.(2.13) does not include a constant term 
with respect to t. Thus we obtain the next proposition, which reduces the amount of cal- 
culation. 

Proposition 2.18. If we choose zeros as integral constants in the formulas (2.11) to 
(2.14) for i = 1, • • • , k - 1, then Eq.(2.13) for i = k is written as 

R k (y) = lim - f G k (s,y, u^(y), ■■■ , uf-%))d S . (2.71) 
f-»°° t J 

If a system for which we calculate the RG equation includes parameters such as k 
and to in Eq. (l2.65l) . then softwares automatically assume that it is in a generic case. 
For example if we compute the RG equation for Eq. (l2.65l) by using Mathematica, it 
is assumed that to ± 0, 1,2. We can verify that the RG equation (12.671) obtained by 
softwares is invalid when to = 0, 1, 2 by computing the second order RG transformation. 
Indeed, it includes the factors to, to - 1, to - 2 in denominators. The explicit form of 
the RG transformation for Eq. fl2.65l) is too complicated to show here. To obtain the RG 
equation for to = 0, 1, 2, substitute them into the original system (12.651) and compute the 
RG equation again for each case. 

3 Restricted RG method 

Suppose that Eq.(2.3) defined on U c M does not satisfy the assumption (B) on U but 
satisfies it on a submanifold No c U. Then, the RG method discussed in the previous 
section is still valid if domains of the RG equation and the RG transformation are re- 
stricted to N . This method to construct an approximate flow on some restricted region is 
called the restricted RG method and it gives extension of the center manifold theory, the 
geometric singular perturbation method and the phase reduction. 

3.1 Main results of the restricted RG method 

Consider the system of the form 

x = f{x) + sg{t, x, s), (3.1) 
defined on U c M. For this system, we suppose that 

(Dl) the vector field / is C°° and it has a compact attracting normally hyperbolic invariant 
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manifold iVo c U. The flow tp t (x) of / on iVo is an almost periodic function with respect 
to t e R uniformly in x e No, the set of whose Fourier exponents has no accumulation 
points. 

(D2) there exists an open set V D No in U such that the vector field g is C 1 in t e R, 
C°° in x e V and small s, and that g is an almost periodic function with respect to t e R 
uniformly in x e V and small s, the set of whose Fourier exponents has no accumulation 
points (i.e. the assumption (A) is satisfied on V). 

If g is independent of t and Eq.(3.1) is autonomous, Fenichel's theorem proves that 
Eq.(3.1) has an attracting invariant manifold N £ near A^ . If g depends on t, we rewrite 
Eq.(3.1)as 

x = f(x) + sg(s, x, s), 

s = 1, K ] 

so that the unperturbed term {fix), 1) has an attracting normally hyperbolic invariant 
manifold R x No on the (s, x) space. Then, a similar argument to the proof of Thm.2.9 
proves that Eq.(3.2) has an attracting invariant manifold iV £ on the (s, x) space, which is 
diffeomorphic to R x A^o- 

In both cases, since N E is attracting, long time behavior of the flow of Eq.(3.1) is well 
described by the flow on iV e . Thus a central issue is to construct _/V e and the flow on it 
approximately. To do so, we establish the restricted RG method on A^o- 

Definition 3.1. For Eq.(3.1), we define C°° maps R h hf : jV -» M to be 

R 1 (y)= lim - [ (p<p s )- l GAs,<pJy))ds, (3.3) 
hl l) (y) = (Dcp t ) y f ((Dcp^G^s, cp s (y)) - R,(y)) ds, (3.4) 

%J oo 



and 



Rfy) = lim - [((Dip^Giis, <p s (y), h™(y), ■■■ , h^(y)) 

1-1 

-(Dtp^YjiDhf^yR^ds, (3.5) 

k=i 

hfiy) = (D(p t ) y f ((Dip^G^tpAyXh^iy),--- ,h { t l) (y)) 

sJ — OO 

i-\ 

-(Dtp^YjiDhf^R^iy) - Ri(y))ds, (3.6) 



i<=i 



for i = 2, 3, • • • , respectively. Note that lim^oo in Eqs.(2.46, 48) are replaced by lim f ^_ 00 
and the indefinite integrals in Eqs.(2.47, 49) are replaced by the definite integrals. With 
these Ri and ny , define the restricted m-th order RG equation for Eq.(3. 1) to be 



y = eRrfy) + s 2 R 2 (y) + ■■■ + s m R m (y), y 6 jV , (3.7) 
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and define the restricted m-th order RG transformation to be 

<*t(y) = <Pt(y) + ^(y) + ■■■ + s m h^(y), y e N , (3.8) 

respectively. 

Note that the domains of them are restricted to N . To see that they make sense, we 
prove the next lemma, which corresponds to Lemma 2.1. 

Lemma 3.2. (i) The maps R t (i = 1,2, • • •) are well-defined and Rfy) e T y N for any 
ye N Q . 

(ii) The maps /z| !) (i = 1, 2, • • • ) are almost periodic functions with respect to t uniformly 
in y e A^o- In particular, Iq are bounded in t e R if y e Nq. 

Proof. Let n s be the projection from T y M to the stable subspace E s of No and n No the 
projection from T y M to the tangent space T y No- Note that n s + n No = id. By the definition 
of an attracting hyperbolic invariant manifold, there exist positive constants C\ and a 
such that 

IWD^vll < d^'IMI, (3.9) 

for any t < 0,y £ Nq and v G T y M. Since (p t (y) e No for all t e R and since G\ is almost 
periodic in t, there exists a positive constant C 2 such that \\Gi(t,(f t (y))\\ < C 2 . Then, 
7T s Ri(y) proves to satisfy 

1 C 

Ik^^y)!! < lim - WnsiDtp^Giis^iyMds 

t->-co t J y 

< lim - f C 1 C 2 e as ds = 0. (3.10) 

t^-oo t J 

This means that R i(y) 6 T y N . 

To prove (ii) of Lemma 3.2, note that the set 

T{S) = {r| WG^s + r^-G^W < 6, \\cp s+T (y)-cp s (y)\\ < S, e R, Vv e iV } (3.11) 
is relatively dense. For t 6 T(§), h\ \ip T (y)) is calculated as 

h { t\y T (y)) = (D<p t \Ay) f (( D <Ps)pl<y)Gi(s, <p s o ^ T (y)) - /?iOp T (y))) 

((IV,.,); 1 ^!^, ^ + r(y)) - (£ty,Wfi(yv(y))) (3.12) 

/ -co 

Putting s' = s + t yields 

Xt+T 
((D<p s ,_ it+T) r y l G 1 (s' - T,<p s ,(y)) - (D(p t ) <PT(y) R 1 ((p T (y)))d S '. (3.13) 
oo 



%J — o 
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Since the space T y N is (D^ f ) v -invariant, n s (D(p t ) y Ri(y) = 0. This and Eq.(3.13) provide 
that 

Xf+T 
\\7r s (Dip sHt+T) )- l \\ ■ ||Gi(5,^(y)) - G^s - T,<p s (y))\\ds 
DO 

Xt+T 
6C ie ais - (,+T)) ds = dCJa. (3.14) 
oo 

Thus we obtain 



where L t , which is bounded in t, is the Lipschitz constant of n s h t |# . This proves that 
tt^ 1 ' is an almost periodic function with respect to t. On the other hand, n NQ h { t l) is written 
as f 

n No h < t l) (y) = n No (D<p t ) y f (^(D^Gi^^Cy))-/?^))^. (3.15) 

U —oo 

Since n No (D<p t ) y is almost periodic, we can show that nM h\ l) (y) is almost periodic in t in 
the same way as the proof of Lemma 2.1 (ii). This proves that h t = n s h t + 7iN h t is 
an almost periodic function in t. The proof of Lemma 3.2 for i = 2, 3, • • • is done in a 
similar way by induction and we omit it here. ■ 

Since Rj(y) 6 T y N , Eq.(3.7) defines a dim/Vo-dimensional differential equation on 

Remark 3.3. Even if h[ are defined by using indefinite integrals as Eqs.(2.47, 49), we 
can show that h[ is bounded as t — > oo, though it is not bounded as t — > -oo. In this case, 
the theorems listed below are true for large t. 

Now we are in a position to state main theorems of the restricted RG method, all 
proofs of which are the same as before and omitted. 

Theorem 3.4. Let a\ m) be the restricted /n-th order RG transformation for Eq.(3.1). If |s| 
is sufficiently small, there exists a function S(t, y, s), y & N such that 

(i) by changing the coordinates as x = a| m) (j), Eq.(3.1) is transformed into the system 

y = eRi(y) + s 2 R 2 (y) + ■■■ + s m R m (y) + s m+l S(t, y, s), (3.16) 

(ii) S is an almost periodic function with respect to t uniformly in y 6 /Vo, 

(iii) S {t, y, s) is C 1 with respect to t and C°° with respect to y e No and s. 

Theorem 3.5 (Error estimate). Let y(t) be a solution of the restricted m-th order RG 
equation and a| m) the restricted m-th order RG transformation. There exist positive con- 
stants so, C and T such that a solution x(t) of Eq.(3.1) with x(0) = a Q m ^(y(0)) e a (No) 
satisfies the inequality 

Mt)-a < r>(y(t))\\<C\s\ m , (3.17) 
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for |e| < £ and < t < T l\e\. 

Theorem 3.6 (Existence of invariant manifolds). Suppose that ^(y) = = Rk-\(y) = 
and e^RjSy) is the first non-zero term in the restricted RG equation for Eq.(3.1). If the 
vector field Rk(y) has a boundaryless compact normally hyperbolic invariant manifold 
L c iVo, then for sufficiently small s > 0, the system (3.2) has an invariant manifold L e on 
the (s, x) space which is diffeomorphic to RxL. In particular, the stability of L £ coincides 
with that of L. 

Theorem 3.7 (Inheritance of symmetries). Suppose that an e-independent Lie group 
H acts on No. If the vector field / and g are invariant under the action of H, then the 
restricted m-th order RG equation for Eq.(3.1) is also invariant under the action of H. 

Recall that our purpose is to construct the invariant manifold N E of Eq.(3.1) and the 
flow on N £ approximately. The flow on N E is well understood by Theorems 3.4 to 3.6, 
and N £ is given by the next theorem. 

Theorem 3.8. Let a* m) be the restricted m-th order RG transformation for Eq.(3 . 1 ) . Then, 
the set {(t, x) \ x € a, (No)} lies within an 0(s m+l ) neighborhood of the attracting invari- 
ant manifold N £ of Eq.(3.2). 

Proof. Though the maps Rfy) and or, are defined on N , we can extend them to the 
maps defined on V D No so that Eq.(3.7) is C 1 close to Eq.(3.16) on V and that No is an 
attracting normally hyperbolic invariant manifold of Eq.(3.7). Then the same argument 
as the proof of Thm.2.9 proves Theorem 3.8. ■ 

If the vector field g is independent of t and Eq.(3.1) is autonomous, we can prove the 
next theorems. 

Theorem 3.9 (Existence of invariant manifolds). Suppose that R\(y) = • • • = 
Rk-i(y) = and e^Rkiy) is the first non-zero term in the restricted RG equation for 
Eq.(3.1) with /-independent g. If the vector field Rtiy) has a boundaryless compact nor- 
mally hyperbolic invariant manifold L, then for sufficiently small s > 0, Eq.(3.1) has 
an invariant manifold L £ , which is diffeomorphic to L. In particular, the stability of L E 
coincides with that of L. 

Theorem 3.10 (Additional symmetry). The restricted RG equation for Eq.(3.1) with t- 
independent g is invariant under the action of the one-parameter group {(p, : _/V — » A^ 1 1 e 
R}. In other words, R t satisfies the equality 

Ri(cp t (y)) = (DtptXRfy), y 6 N , (3.18) 

for i = 1, 2, • • • . 

For autonomous systems, Thm.3.8 is restated as follows: Recall that if the function 
g depends on t, the attracting invariant manifold Af £ of Eq.(3.1) and the approximate 
invariant manifold described in Thm.3.8 depend on t in the sense that they lie on the 
(s, x) space. If Eq.(3. 1) is autonomous, its attracting invariant manifold Af £ lies on M and 
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is independent of t. Thus we want to construct an approximate invariant manifold of N E 
so that it is also independent of t. 

Theorem 3.11. Let ar| m) be the restricted m-th order RG transformation for Eq.(3.1). If g 
is independent of t, the set a] (No) = {of 1 ' (y) \ y e N } is independent of t and lies within 
an 0(s m+l ) neighborhood of the attracting invariant manifold N £ of Eq.(3.1). 

Proof. We have to show that the set a| m) (/V ) is independent of t. Indeed, we know the 
equality h { t '\iff(y)) = h® f (y) as is shown in the proof of the Thm.2.15. This proves that 

aJJCiVb) = a^\MNo)) = ^(No). (3.19) 
The rest of the proof is the same as the proofs of Thm.2.14 and Thm.3.8. ■ 

3.2 Center manifold reduction 

The restricted RG method recovers the approximation theory of center manifolds (Carr 
[7]). Consider a system of the form 

x = Fx + sg(x, s) 

= Fx + egi(x) + s 2 g 2 (x) + ■■■ , xe R", (3.20) 

where unperturbed term Fx is linear. For this system, we suppose that 

(El) all eigenvalues of the n x n constant matrix F are on the imaginary axis or the left 
half plane. The Jordan block corresponding to eigenvalues on the imaginary axis is diag- 
onalizable. 

(E2) g is C°° with respect to x and s such that g(0, e) = 0. 

If all eigenvalues of F are on the left half plane, the origin is a stable fixed point and the 
flow near the origin is trivial. In what follows, we suppose that at least one eigenvalue 
is on the imaginary axis. In this case, Eq.(3.20) has a center manifold which is tangent 
to the center subspace /V at the origin. The center subspace N Q , which is spanned by 
eigenvectors associated with eigenvalues on the imaginary axis, is an attracting normally 
hyperbolic invariant manifold of the unperturbed term Fx, and the flow of Fx on N is 
almost periodic. However, since N is not compact, we take an n-dimensional closed 
ball K including the origin and consider iVo n K. Then, we obtain the next theorem as a 
corollary of Thm.3.9 and Thm.3. 1 1 . 

Theorem 3.12 (Approximation of Center Manifolds, [12]). Let or, be the restricted 
m-th order RG transformation for Eq.(3.20) and K a small compact neighborhood of the 
origin. Then, the set a) (KDNq) lies within an 0(e" I+1 ) neighborhood of the center man- 
ifold of Eq.(3.20). The flow of Eq.(3.20) on the center manifold is well approximated by 
those of the restricted RG equation. In particular, suppose that Ri(y) = ■ ■ ■ = Rk-i(y) = 
and e k Rk(y) is the first non-zero term in the restricted RG equation. If the vector field 
Rk(y) has a boundaryless compact normally hyperbolic invariant manifold L, then for 
sufficiently small e > 0, Eq.(3.20) has an invariant manifold L s on the center manifold, 
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which is diffeomorphic to L. The stability of L E coincides with that of L. 
See Chiba [12] for the detail of the proof and examples. 



3.3 Geometric singular perturbation method 

The restricted RG method can also recover the geometric singular perturbation method 
proposed by Fenichel [19]. 

Consider the autonomous system 

x = f(x) + sg(x, s), xeR" (3.21) 

on R" with the assumption that 

(F) suppose that / and g are C°° with respect to x and s, and that / has an /^-dimensional 
attracting normally hyperbolic invariant manifold N Q which consists of fixed points of /, 
where m < n. 

Note that this system satisfies the assumptions (Dl) and (D2), so that Thm.3.9 to 
Thm.3.11 hold. The invariant manifold N consisting of fixed points of the unperturbed 
system is called the critical manifold (see [1]). For this system, Fenichel[19] proved that 
there exist local coordinates (u, v) such that the system (3.21) is expressed as 

( u = sg x (u, v, s), u G R m , 

\ v = f(u, v) + sg 2 (u, v, s), v e R"~ m , { ' ; 

where f(u, 0) = for any u 6 R m . In this coordinate, the critical manifold A^o is locally 
given as the w-plane. Further he proved the next theorem. 



Theorem 3.13 (Fenichel [19]). Suppose that the system u = sgi(u,0, 0) has a compact 
normally hyperbolic invariant manifold L. If e > is sufficiently small, the system (3.22) 
has an invariant manifold L E which is diffeomorphic to L. 

This method to obtain an invariant manifold of (3.22) is called the geometric singular 
perturbation method. By using the fact that (p t (u) = u for u e A^ , it is easy to verify 
that the system u = sg\{u, 0, 0) described above is just the restricted first order RG equa- 
tion for Eq.(3.22). Thus Thm.3.13 immediately follows from Thm.3.9. Note that in our 
method, we need not change the coordinates so that Eq.(3.21) is transformed into the 
form of Eq.(3.22). 

Example 3.14. ©Consider the system on R 2 

| X\ = -X\ + (xj + c)x2, 
\ sx 2 = x\ - (xi + Y)x z , 
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where < c < 1 is a constant. This system arises from a model of the kinetics of enzyme 
reactions (see Carr [7]). Set t = ss and denote differentiation with respect to s by ' . 
Then, the above system is rewritten as 

Xj = s(—X\ + CX2 + X1X2), ^ 

x'2 = X\ — X2 ~ X\X2- 

The attracting critical manifold No of this system is expressed as the graph of the function 

x 2 = Hxi) ■= T~ • (3-25) 

1 + X\ 

Since the restricted first order RG transformation for Eq.(3.24) is given by 

a! '' fri Hw<r+JH-< c -iw<> <3 - 26) 

Theorem 3.11 proves that the attracting invariant manifold of Eq.(3.24) is given as the 
graph of 

X\ (c — l)jCi 7 

\Tx~ l ~ S {\ +Xi) 4 
If \x\ I is sufficiently small, it is expanded as 

x 2 = X\{\ - X]) - s{c - \)x\{\ - Ax\) + 0(x\, s 2 ) 

= (1 -s(c- l))xi -(1 -4s(c- l))x 2 + 0(x\,s 2 ). (3.28) 

This result coincides with the result obtained by the local center manifold theory (see 
Carr [7]). The restricted first order RG equation on /Vo is given by 

y x = s—— . (3.29) 

1 + y x 

This RG equation describes a motion on the invariant manifold (3.27) approximately. 
Since it has the stable fixed point y\ = if c < 1, the system (3.24) also has a stable fixed 
point (xi,x 2 ) = (0, 0) by virtue of Theorem 3.9. 



x 2 = - — — - e \ , - ^ + 0{e 2 ). (3.27) 



3.4 Phase reduction 

Consider a system of the form 

x = f(t, x) + £gi(t, x), x 6 R". (3.30) 
For this system, we suppose that 

(Gl) the vector fields / and g are C°° in x, C 1 in t, and T-periodic in t. It need not be the 
least period. In particular, / and g are allowed to be independent of t. 
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(G2) the unperturbed system x = f{t, x) has a fc-parameter family of r -periodic solutions 
which constructs an attracting invariant torus T* c R". 

Let a = (ari,--- ,at) be coordinates on T k , so-called phase variables. It is called 
the phase reduction to derive equations on a which govern the dynamics of Eq.(3.30) 
on the invariant torus. The phase reduction was first introduced by Malkin [33,34] and 
rediscovered by many authors. In this subsection, we show that the RG method can 
recover the phase reduction method. 

The next theorem is due to Malkin [33,34]. See also Hahn [23], Blekhman [4], and 
Hoppensteadt and Izhikevich [25]. 

Theorem 3.15 (Malkin [33,34]). Consider the system (3.30) with the assumptions (Gl) 
and (G2). Let a = (a\, ■ ■ • , at) be phase variables and U(t; a) the periodic solutions of 
the unperturbed system parameterized by a. Suppose that the adjoint equation 

--(^f(t, U(t;a))\ Qi (3.31) 



dt \dx 

has exactly k independent 7 -periodic solutions Qi(t; a), ■ ■ ■ , Qk(t; a), where A T denotes 
the transpose matrix of a matrix A. Let Q = Q(t; a) be the k x n matrix whose columns 
are these solutions such that 

Q J —(t;a) = id. (3.32) 
oa 

Then, Eq.(3.30) has a solution of the form 

x(t) = U(t, a{t)) + 0(s), (3.33) 
where a(t) is a solution of the system 

^ = ^f Q(s;a) T gl ( S ,U(s;a))d S . (3.34) 
dt T Jo 

Now we show that the system (3.34) of the phase variables is just the first order 
RG equation. Note that the system (3.30) satisfies the assumptions (Dl) and (D2) with 
A^o = and the RG method is applicable. The restricted first order RG equation for 
Eq.(3.30) is given by 

y = s \i m f(^O0) gl ( S ,cp s (y))ds, yeT k . (3.35) 
Jo \ oy ) 

Let us change the coordinates by using the ^-parameter family of periodic solutions as 
y = U(0; a). Then, Eq.(3.35) is rewritten as 

dU /n x . f*(d<p t 



— (0;a)a = elim Up(tf(0;a)) gl (s,U(s;a))ds. (3.36) 
da f->-« Jo \ dy 
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Since U(t; a) = (p t (U(0; a)), the equality 




(3.37) 



holds. Then, Eqs.(3.32), (3.36) and (3.37) are put together to obtain 



a = s lim 



= s lim 



= s lim 



?— >— oo 




(3.38) 



Since the integrand in the above equation is ^-periodic, Eq.(3.38) is reduced to Eq.(3.34). 

4 Relation to other singular perturbation methods 

In the previous section, we have seen that the restricted RG method unifies the center 
manifold reduction, the geometric singular perturbation method and the phase reduction. 
In this section, we show that the RG method described in Section 2 unifies the tradi- 
tional singular perturbation methods, such as the averaging method, the multiple time 
scale method and the normal forms theory. We will also give explicit formulas for reg- 
ular perturbation solutions and show how the RG equation is derived from the regular 
perturbation method through the envelope theory. 

4.1 Averaging method 

The averaging method is one of the most traditional and famous singular perturbation 
methods based on an idea of averaging a given equation, which is periodic in time t, 
with respect to t to obtain an autonomous equation (see Sanders, Verhulst and Murdock 
[47]). In most literature, only first and second order averaging equations are given and 
they coincide with first and second order RG equations, respectively. Our m-th order RG 
equation gives a generalization of such lower order averaging equations. In Chiba and 
Paz6[13], the third and fifth order RG equations are used to unfold a phase diagram of 
the Kuramoto model of coupled oscillators. 

Many authors formulated the averaging method for time-periodic differential equa- 
tions. However, the RG equation (the averaging equation) can be defined as long as 
limits in Eqs.(2.11, 13) exist and Rfy) are well-defined even if a system (2.1) does not 
satisfy the assumption (A). If Ri, ■ ■ ■ ,R m are well-defined for a system (2.1), which need 
not satisfy the assumption (A), we say (2.1) satisfies the KBM condition up to order m 
(KBM stands for Krylov, Bogoliubov and Mitropolsky, see Bogoliubov and Mitropolsky 
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[6]). We showed that a system (2.1) with the assumption (A) satisfies the KBM condition 
up to all order (Lemma 2.1). If the assumption (A) is violated, the error estimate (2.27) 
for approximate solutions may get worse, such as \\x(t) - a| m) (y(0)ll < C ^je, even if the 
KBM condition is satisfied and thus the RG equation is defined. See Sanders et al. [47] 
and DeVille et al. [32] for such examples. 

Even if (2.1) satisfies the KBM condition, we need additional assumptions to prove 
Thm.2.9. It is because to apply Fenichel's theorem, which was used in the proof of 
Thm.2.9, we have to show that the error function S (t, y, s) in Eq.(2.23) is C 1 with respect 
to t, x and bounded in t e R. See Chiba [10] for more detail. 



4.2 Multiple time scale method 



The multiple time scale method [3,38,44] was perhaps introduced in the early 20th cen- 
tury and now it is one of the most used perturbation techniques along with the averaging 
method. In this subsection, we give a brief review of the method and show that it yields 
the same results as the RG method. 

Consider the system (2.1) satisfying the assumption (A). Let us introduce the new 
time scales t m as 

t m = s m t, m = 0,1,2,- ■■ , (4.1) 

and consider t and to, h,t 2 ,-- • to be independent of each other. Then, the "total" deriva- 
tive d/dt is rewritten as 



d 



d 



d 



+ S— + £ + • • • . 

dt dt dt\ dt 2 
Let us expand the dependent variable x as 

2 

X — Xq + £X\ + S X 2 + " " " , 

where x t = x,(? , h, t 2 , • • ■ ). Substituting Eqs.(4.2) and (4.3) into (2.1), we obtain 



(4.2) 
(4.3) 



(ir + € ^r + e2 ^r + ) ( Xo + €Xi + £2xi + ) = ^ ^s^o, *o + + s 2 x 2 +•••)• 

(4.4) 

Expanding the both sides of the above in e and equating the coefficients of each e k , we 
obtain ODEs on x , xi, x%, ■ ■ ■ ; 



= 0, 



8xq 
dt 

dx y dx 

—— — (J\\tQ, Xq) — 



dt 
dx 2 



dh 



— G 2 (tg, Xq, Xi) 



dxi dxo 
~dh~~dT 2 



(4.5) 



— G m (t , Xq, ■ ■ ■ , X OT -l) — 



dx ■ 
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where the functions G^, k = 1,2, ••• are defined through Eq.(2.6). Let x = y = 
y(h,h, ■ ■ ■ ) be a solution of the zeroth order equation. Then, a general solution of the 
first order equation of (4.5) is given as 

x 1 =B 1 + J (g^v) - ^]ds = B 1 + u%(y) + R^to - $Lto, (4.6) 

where B\ = B[(y; t\, t 2 , ■ ■ ■ ) is independent of to and where Ri and are defined by 
Eqs.(2.11) and (2.12), respectively. Now we define dyjdt\ so that x\ above is bounded in 

to : 

Ri(y)t - % = 0. (4.7) 

Ot\ 

This condition is called the non-secularity condition and it yields 

^ = R x (y), x x = Xl (y) = B l+ (4.8) 
Next thing to do is calculating x 2 . The equation on x 2 is written as 

f* = G2(f0 ,,, Bl+B «V))-| : (B 1+ <: , w)-| 

dgi m du ?J dB l dB x dy 

= ■^■(fo.yX^i + «<>)) + g 2 (t ,y) - -^-(yWi(y) - - -^iO0 - ^4.9) 

a general solution of which is given by 

f°(dgi m du { p dBi dBx dy\ 

x 2 = B 2 + J I^^X^+^^ + ^^-^^^^-^-^Kv)-^!^ 



cfy <9fi dy dt2 



n\ du tl 
B 2 + u%(y) + ~^(y)Bi 



+-±(y)B x to + R 2 (y)t Q - ^t Q - —±Ry(y)to - ^t Q , (4.10) 
oy ot\ oy ot 2 

where B 2 = B 2 (y; t\,t 2 , ■ • • ) is independent of to. If we impose the non-secularity condi- 
tion so that x 2 is bounded in to, we obtain 

dy dBi 

= Ri(y) t~ — [B u Rim, (4.11) 

Ot 2 Ot\ 

du (l) 

x 2 = x 2 (y) = B 2 + «®(y) + -^-Pi, (4.12) 
where the commutator [•, •] is defined in Eq. (l2.61l) . 
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By proceeding in a similar manner, we can show that the non-resonance condition at 
each stage yields equations on y of the form 

% = Rfy) + Qi(Ru ■ ■ ■ ,Ri-i,B u • • • , Bi- 2 )(y) - ^ - [^-i^ilCv), (4.13) 

Oti Ot\ 

for i = 2, 3, • • • , where Qj is a function of R\, • • • and B\, ■ • ■ , B^ 2 . Combining these 
equations on y, we obtain 

dy dy dy 2 dy 

= — h £— 1- S 1 

dt dt ot\ dt 2 

= + s 2 ^R 2 (y)- ^ - [B u R{\(yjj + ••• . (4.14) 

The functions B\,B 2 , - ■ • are still left undetermined. We propose two different ways to 
determine them. 

(i) If we define 5,'s as solutions of the linear inhomogeneous differential equations 

^ = -[B u R{\<y) + R M <y) + Q M (Ru ■ ■ ■ ,R h B u ■ ■ ■ ,fi M )(y), (4-15) 
then Eq.(4.14) is reduced to the equation 

^=sR 1 (y). (4.16) 
dt 

Let y = y(t) be a solution of this equation. Then, we can prove that a curve defined by 

it)) + ■■■ + s m x m (y(t) 
( a„(D 



at n) (y(t)) := y(t) + £ x 1 (y(t)) + £ 2 x 2 (y(t)) + --- + £ m x m (y(t)) 



= y + s{B x + i^OO) + s 2 



B 2 + u (2 \y) + ^-(y)B l 
dy 



+ 



(4.17) 

y=y(i) 



provides an approximate solution for Eq.(2.1) and satisfies the same inequality as (2.27) 
(see Murdock[38]). 

(ii) Otherwise, we suppose that 5,'s are independent of h,t 2 ,--- so that Eq.(4.14) is 
independent of t\, t 2 , ■ ■ ■ . Then Eq.(4.14) takes the form 

% = sRtiy) + s 2 (R 2 (y) - [B^R^y)) + • • • , (4.18) 
dt 

which coincides with the (simplified) RG equation (see Sec. 2. 4). 
4.3 Normal forms 

Let us consider the system (2.5) with assumptions (CI) and (C2). The technique of 
normal forms is used to analyze local dynamics of such a system and had been well 
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developed in the last century (Murdock [39]). Let us recall the definition of normal forms. 
For a given system (2.5), there exist a time-independent local coordinate transformation 
x = h(z) defined near the origin which brings Eq.(2.5) into the form 

z = Fz + sg x (z) + ■■■ + s m ~g m {z) + s m+l S(z,s) (4.19) 

with the properties that gj(z), i= 1, • • • , m satisfy gj(e Ft z) = e F 'gj(z) and that S is C°° in z 
and e. Then the truncated system 

z = Fz + egi(z) + • • • + s m g m (z) (4.20) 

is called the normal form ofEq.(2.5) up to order m. 

Note that if the matrix F is not diagonalizable, different definitions are adopted for 
normal forms. 

Now we consider the RG equation and the RG transformation for Eq.(2.5), respec- 
tively, defined by Eqs.(2.50) and (2.51) with Eqs.(2.46) to (2.49). The flow <p t (y) in 
Eqs.(2.46) to (2.49) is given as (p t (y) = e Ft y by using the fundamental matrix e Ft in our 
situation. Thm.2.5 shows that the RG transformation x = a^ m) (y) brings Eq.(2.5) into 
the form of Eq.(2.23). Because of Thm.2.15, further change of variables as y = e~ Ft z 
transforms Eq.(2.23) into the system 

z = Fz + sRriz) + ■■■ + e m R m {z) + e m+l e F, S (t, e~ Ft z, s) (4.21) 

with the property that Rj(e Ft z) = e Ft Ri{z) for i = 1, • • • , m. 
Now we need the next lemma. 

Lemma 4.1. For the autonomous system (2.5), the RG transformation satisfies the equal- 
ity a ( ; n \e Ft 'y) = a™(y) for all t, f,y. 

This lemma immediately follows from the equality h^\(p f {y)) = hf +t ,{y), i = 1, 2, • • • 
proved in the proof of Thm.2.15. Thus it turns out that the composition of two transfor- 
mations x = 0Cf(y) and y = e~ F, z is independent of t: 

x = a ( t m \e- Ft z) = af\z). (4.22) 

This proves that Eq.(4.21) is obtained from Eq.(2.5) by the time-independent transforma- 
tion, and thus the truncated system 

z = Fz + sRi (z) + ■ ■ ■ + s' n R m (z) (4.23) 



is a normal form of Eq.(2.5) up to order m. 

Note that normal forms of a given system are not unique in general and the simplest 
form among them is called a hyper-normal form ([2,39,40]). The RG method can also 
provide hyper-normal forms by choosing undetermined integral constants in Eqs.(2.47) 
and (2.49) appropriately as was discussed in Sec. 2.4. See also Chiba [11]. 
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4.4 Regular perturbation and envelopes 



In this subsection, we give some properties of the regular perturbation method. The RG 
equation is obtained from regular perturbation solutions through Kunihiro's idea based 
on envelopes [28,29]. 

Consider the system (2.1) on R" with the assumption (A). Let us construct a formal 
solution of it of the form 



x 



Xq + SX\ + S X2 + 



(4.24) 



Substituting Eq.(4.24) into Eq.(2.1) and equating the coefficients of each s k , we obtain a 
system of ODEs : 

( x = 0, 

X\ = Gi(t, Xq), 



(4.25) 



Xk — Gk(t, Xq,- - • , Xk-\), 



where the functions Gk are defined through Eq.(2.6). Solving these equations and substi- 
tuting solutions Xk = xt(t) into Eq.(4.24), we obtain a formal solution x = x(t). If Eq.(2.1) 
is analytic in e e £ c C, Eq.(4.24) converges for s e £ and gives an exact solution of 
Eq.(2.1). However, in this section, we need not such an assumption and regard Eq.(4.24) 
as a formal power series in s. To construct a formal solution x(t) as above is called the 
regular perturbation method. 

Let us define Rfy) and u [ t l) (y), i = 1,2, • • • by Eqs.(2.11) to (2.14). By using these 
functions, regular perturbation solutions (4.24) are given as follows: 

Proposition 4.2. Solutions of the system of ODEs (4.25) are given as 



x k = x k (t,y) = uf(y) + p™(t,y)t + p\ K, (t,y)t z + ■■■+ pf(t,y)f, 



ft) 



(10, 



,(*)/ 



(4.26) 



for k = 1, 2, • • • , where y 6 R" is an arbitrary constant (a solution of xq = 0) and where 
pY> are defined by 



p\ 1 \t,y) = p\ L, (y)=Ri(y) 



pf(t,y) 



pf(t,y) 



pf(t,y) = p)"(y) = - 



duf 



J = 2,3,. 



J 



k=i 
9P% 



k=\ 



dy 



(4.27) 



(0/ 



l dpi, 



{ P f(t,y) = 0, 



i dy 



(y)Ri(y), i = 2,3, 
j > i. 



In particular, py(t,y) are almost periodic functions with respect to t for all i, j e N. 
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Proof. This proposition is proved in Chiba [10]. 



Now we derive the RG equation (2.18) according to an idea of Kunihiro [28,29]. We 
can show that regular perturbation solutions are also expressed as 



x = x(t, r,y) = y + sxi(t, r,y) + s x 2 (t, r, y) + 



(4.28) 



with 



x k (t, t, y) = ufiy) + p\ K, (t, y)(t - r) + p%>(t, y){t — t) l + • • • + pf(t, y){t - r)\ (4.29) 

by choosing initial times and initial values appropriately when solving the system (4.25) 
so that a formal solution x(t, r,y) passes through y + ££Lj €^uf\y) att = T, where r e R 
is an arbitrary constant. Further, we regard y = y(r) as a function of r to be determined, 
and consider the family {3c(t, r, y(r))} Te R of regular perturbation solutions parameterized 
by r e R. It is known that regular perturbation solutions are close to exact solutions only 
for a short time interval \t\ ~ 0(1). However, if we move the initial value y + 2£i s*Mf'(y) 
along an exact solution by varying r as Fig .7, it seems that the envelope of the family 
{x(t,T,y(T))} re R gives the exact solution of Eq.(2.1). 



Ak) 



(*)/ 



exact solution 




regular perturbation solutions 



Fig. 7: Exact solution of Eq.(2.1) and a family of regular perturbation solutions 
x(t,T,y(T)). 

The envelope is calculated as follows: At first, we differentiate x(t,T,y(r)) with re- 
spect to t = t and equate it to zero : 

d 



dr 



x(t,T,y(T)) = 0. 



Substituting Eqs.(4.28) and (4.29) into the above yields 



rlt £—1 rlr 



dt dr 



x k (t,T,y(r)) 



duf\ dy 



k=i 



k=\ 



(4.30) 



(4.31) 



41 



This equality is formally satisfied if we define y = y(t) as a solution of the ODE 

, CO 

-£=yy**oo- (4.32) 

This recovers the m-th order RG equation if truncated at order s m . With this y(t), the 
envelope for the family {x(t, r,y(r))} TeR is given by y(t)). Again Eqs.(4.28) and 
(4.29) prove that 

oo 

x(t,t,y(t)) = yiO + J^^x^^yit)) 
k=i 

oo 

k=\ 

= a ( r\y(t)) + 0(s m+ '), (4.33) 
which gives the m-th order RG transformation if truncated at order e m . 

In the rest of this section, we prove that a formal solution (4.24) with Eq.(4.26) ob- 
tained by the regular perturbation method includes infinitely many convergent subseries 
even if Eq.(2.1) (and thus Eq.(4.24)) is not analytic in e. 

Proposition 4.3. For sufficiently small \st\, the series £/=i ^t l p] (t,y) are convergent 
farJfc = 0,l,2,-.-. 

To prove Prop.4.3, we need the next simple lemma. 

Lemma 4.4. Let h\(t) and h 2 (t) be almost periodic functions with Mod^) d Mod(/? 2 ) 
(the module of almost periodic functions is defined in Sec. 2.1). If 

Km(h l (t)-h 2 (t)) = 0, (4.34) 

t— »co 

then hi(t) = h 2 (t) for all t e R. 

Proof of Lemma 4.4. Functions hi, i = 1,2 are expanded in Fourier series as hi(f) = 
E*eM*K*i) Ci(A k )e iAkt . Then, Eq.(4.34) is rewritten as 

lim V ( Cl (^)-c 2 (4))^ = 0. (4.35) 

A t eMod(fti) 

This proves that C\(Au) = c 2 (/LJ for all A k e Mod^x). ■ 
Proof of Prop.4.3. Let 



x(t) = y Q + f j eWpfiy,) + £ A uf(y Q ) + £ sVrf^ y ) 



k=l v /=1 



(4.36) 
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be a formal solution of Eq.(2.1) constructed by the regular perturbation method and let 
y = y(st, s) be a solution of the m-th order RG equation 



dy 
d(et) 



= R l (y) + sR 2 (y) + --- + e m ~ l R m (y) 



(4.37) 



for Eq.(2.1). By Thm.2.7, there exist positive constants C and T such that the inequality 
\\x(t) - a ( r\y(st, e))\\ < C\s\ m , 0<t<T/s (4.38) 
holds if we choose initial values of x(t) and y(st, s) appropriately. Putting t = T/s yields 

m 

af/ s (y(T, e)) = y(T, s) + J] ^4] s (y(T, e)) 

k-i 

co m—1 / oo \ 

= ?o + ^VrfW + £V uf^ + ^T'pf^iT/e^o) \ + 0(s m ). (4.39) 



Note that functions uf\y) and p y j +K> (t,y) are almost periodic functions in t such that 



k=\ v 



a+k) 



i=i 



Mod(wf ) ) and Mod(p)" KI ) are included in Mod(g). In particular, they are bounded in 



a+k) 



t e R. Thus taking the limit e — » in Eq.(4.39) provides 

oo 

y(.T,0)=y Q + Y j T l p ( I \y ). 



(4.40) 



/=i 



This proves that the series T l pf(yQ) is convergent. 

Next thing to do is to prove Prop.4.3 for k = 1 . Using Eq.(4.40) and dividing the both 
sides of Eq.(4.39) by s, we obtain 



m—\ 



kl 

Taking the limit e — » yields 



k=\ 



i=i 



+0{s m ~ l ). 



dy 

OS e->0 



u%(y(T,0)) - 4l(y ) - J] T l pf +1) (T /s,y ) 
< i=i 

Now Lemma 4.4 proves that 

oo „ 

J] T l pf +1) (t,y ) = f-(T,0) + ^ !) (y(r,0)) - u^(y ), 



0. 



(4.41) 



(4.42) 



(4.43) 



i=i 

and the series J^i T l p\ (t,yo) proves to be convergent. 

The proof of Prop.4.3 for k > 2 is done in a similar way and omitted. 
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Though it seems that Prop. 4. 3 has no relationship to the RG method, it is instructive to 
notice the equality (4.40). Since y(T, 0) in the left hand side is obtained from a solution 
of the first order RG equation dy/d(st) = Ri(y), the series 2/=i T l p^(yo) proves to be 
determined by only the first order RG equation. This is remarkable because if we want to 
obtain T l p^(yo) by using the regular perturbation method, we have to solve infinitely 
many ODEs (4,25). A similar argument shows that a solution of the m-th order RG 
equation involves all terms of the form e^ +1 t'p ( l +k \t,y Q ) fork = 0, • • • ,m - 1 and / = 
1, 2, • • • in the formal solution (4.24). 



5 Infinite order RG equation 

Infinite order RG equations and transformations are not convergent series in general. In 
this section, we give a necessary and sufficient condition for the convergence. It will be 
proved in Sec. 5. 3 that infinite order RG equations for linear systems are convergent and 
related to Floquet theory. 

5.1 Convergence condition 

Let us consider the system 

x = sg(t, x, s), xgM (5.1) 
on a real analytic manifold M with the following assumption : 

(A') The vector field g is analytic with respect to t e R, x e M and s £ I, where / c R is 
an open interval containing 0. Further g is T -periodic in t. 

Because of the periodicity, we can regard (5.1) as a system on S 1 x M : 

-l s \-( 1 

dt \ x ) ~ \ sg(s, x, s) 

where 5 1 is a circle with a C w structure. 

Theorem 2.5 states that the infinite order RG transformation defined by 

x = a t (y) :=y + su ( t l \y) + s 2 uf ] (y) + ■■■ (5.3) 

formally brings the system (5.1) into the infinite order RG equation 

y = sR l (y) + s 2 R 2 (y) + -- - , (5.4) 

where "formally" means that Eqs.(5.3) and (5.4) are not convergent in general. If Eq.(5.3) 
is convergent, so is Eq.(5.4). A necessary and sufficient condition for the convergence of 
Eq.(5.3) is given as follows: 

Theorem 5.1. For the system (5.1) with the assumption (A'), there exist an open neigh- 
borhood U = U y of S 1 x {y} x {0} in S 1 x M x / for each y e M and an analytic infinite 



(5.2) 
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order RG transformation on U, if and only if the system (5.1) is invariant under the T 1 
(1 -torus) action of the form 

T 1 : (t, x) i-» (t + k, x + scr k (t, x, e)), k e R, (5.5) 

where cr k (t, x, s) is analytic with respect to k, t, x, s and T-periodic in k and t. 

Recall that RG equations and RG transformations are not unique and not all of them 
are convergent even if the condition of Thm.5.1 is satisfied. 

The proof of this theorem involves Lie group theory and will be given in Sec. 5. 2, 
though the idea of the proof is shown below. 

Since the infinite order RG equation is an autonomous system, it is invariant under 
the translation of t, (t, y) i-» (t + k, y). If an infinite order RG equation and transformation 
are convergent and well-defined, the system (5.1) is invariant under the action defined by 
pulling back the translation by the RG transformation : 

(t,x)h* (t + k,a t+k o a ; l (x)). (5.6) 

Since a, is r -periodic (Lemma 2.1 (ii)), this defines the T 1 action on the space S 1 x M. 
Conversely, if the system (5.1) is invariant under the action (5.5), then a simple extension 
of Bochner's linearization theorem proves that there exists a coordinates transforma- 
tion x i — > y such that the action (5.5) is written as (t, y) h-» (t + k, y). We can show that this 
transformation is just an RG transformation. See Section 5.2 for the detail. 

In the rest of this subsection, we consider an autonomous system on C" of the form 

x = Fx + sg(x, s) 

= Fx + s gl (x) + s 2 g 2 (x) + ■■■ , xeC", (5.7) 

with the assumptions (C1),(C2) (see Sec.2.1) and (C3) below : 
(C3) g(x, s) is analytic with respect to x e C" and e e / c R. 

For this system, RG transformations are defined by Eq.(2.51) with Eqs.(2.46) to 
(2.49). The next corollary immediately follows from Thm.5.1. 

Corollary 5.2. Suppose that all eigenvalues of F have pairwise rational ratios. Then, 
there exist an open neighborhood U = U y of S 1 x {y} x {0} in S 1 x M x / for each y e M 
and an analytic infinite order RG transformation on U, if and only if the system (5.7) is 
invariant under the T 1 action of the form 

T 1 : x h-> e Fk x + scr k {t, x, s), k 6 R, (5.8) 

where cr k (t, x, s) is analytic with respect to k, t, x, e and periodic in k and t. 

Proof. By changing the coordinates as x = e Ft X, Eq.(5.7) is rewritten as X = ee~ Ft g(e Ft X, s). 
Since e Ft is periodic because of the assumption of Corollary 5.2, we can apply Thm.5.1 
to this system. Note that we do not need the assumption (C2) to prove Corollary 5.2. ■ 
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Recall that RG equations for Eq.(5.7) are equivalent to normal forms (Sec. 4. 3). If 
there are irrational ratios among eigenvalues of F, Thm.5.1 is no longer applicable. For 
such a system, Zung [54] gives a necessary and sufficient condition for the convergence 
of normal forms of infinite order, although he supposes that gi(x) in Eq.(5.7) is a homo- 
geneous vector field of degree i for i = 2, 3, • • • . A necessary and sufficient condition for 
the convergence of infinite RG equations and transformations for Eq.(5.7) is given in a 
similar way to Zung's theorem as follows: 

Remember that an RG equation for Eq.(5.7) has the property that Rj(e Fk y) = e Fk Ri(y), k e 
R for i = 1, 2, • • • , if integral constants in Eqs.(2.47) and (2.49) are appropriately chosen 
(Thm.2.15). If a matrix B satisfies the equalities FB = BF and Q(e Bk y) = e Bk Q(y), k e R 
for all polynomial vector fields Q such that Q(e Fk y) = e Fk Q(y), k e R, then B is called 
subordinate to F. It is known that if F is a diagonal matrix, we can take B having the 
form B = diag {ib\, • • • , ib n ), where i = V-T and bj e Z for j = 1, • • • ,n (see Murdock 
[39], Zung [54]). Let p be the maximum number of linearly independent such matrices 
B { , - ■■ ,B p and call it the toric degree of F. Then, the matrix e B < k i+-+ B P k i^ (fc lf ■ • ■ , k p ) e 
R p induces the T p action on C" and the RG equation is invariant under this action. By 
a similar way to the proof of Thm.5.1, we can prove the next theorem, whose proof is 
omitted here. 

Theorem 5.3. Let p be the toric degree of F. For the system (5.7), there exist an open 
neighborhood U = U y of S 1 x {y} x {0} in S 1 x M x / for each y e M and an analytic 
infinite order RG transformation on U, if and only if (5.7) is invariant under the T p action 
of the form 

TP-.x^ e Blkl+ - +B ^x + eo- kl ,., kp (t, x, e), (h, ■ ■ ■ , k p ) e R p , (5.9) 

where cr^ ... kp (t, x, s) is analytic with respect to k\ , • • • ,k p , t, x, e and periodic in k\ , • • • ,k p 
and t. 

5.2 Proof of Theorem 5.1 

In this subsection, we give a proof of Thm.5.1. At first, we provide a few notations and 
facts from Lie group theory. 

Let K be a compact Lie group, dk a bi-invariant measure on K, V a complete locally 
convex topological vector space, and n a representation of K in V. We define the average 
av(^) : V — » V of n by 

av(7r)(v) = f n(k)(v)dk. (5.10) 

JK 

Then, the next theorem holds (see Duistermaat and Kolk [15] for the proof). 

Theorem 5.4. The average av(7r) is a linear projection from V onto the space : = 
{v 6 V | n(k)v = v, Vk e K}; that is, equalities 

n(k) o av(7r)(v) = av(7r)(v) (5.11) 
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and 

av(7r)(v) = v, veV n(K) (5.12) 

hold. 

Let S 1 be a circle with a C w structure, M a real analytic manifold, and /cRan open 
interval containing 0, as Sec.5.1. The next theorem is a simple extension of Bochner's 
linearization theorem [5,15]. 

Theorem 5.5. Let A : K -» DifP(S 1 x M x /) be a C action of K on S 1 x M x / such 
that A(k) is expressed as 

A(k)(t, x, s) = (ai(k)(f), a 2 (k)(t, x, s), e), k 6 K, (5.13) 

where ai(fc) : S 1 — » 5 1 and a 2 (£) : S 1 x M x / — > M are C w maps. Suppose that there 
exists xq e M such that 

a 2 (k)(t,xo,0) = x . (5.14) 

Then, there exist an open neighborhood U of S 1 x {jc } x {0} in 5 1 x M x / and a C w 
diffeomorphism ^ from U into 5 1 x T XQ M x R such that 

Tp o A(k) = (ai(k) x D x a 2 (k)(t, x , 0) x id) o ^, (5.15) 

where Z) v is the derivative with respect to x. The <£> is expressed as 

Tp{t, x, s) = (t, a-sit, x, e), s) (5.16) 

with a C w map a 3 and satisfies 

£(f, jc , 0) = (t, 0, 0), DJpit, x Q , 0) = id. (5.17) 

Further if we suppose 

D n x a 2 (k)(t,x Q ,0) = (n>2), (5.18) 

then satisfies 

D^(?,x ,0) = (n>2). (5.19) 

Proof. Let W be a .fif- invariant open neighborhood of S 1 x {xo} x {0} and V a space of C w 
maps v : ^^S'xJ ri MxR such that 

v(5 1 x {x } x {0}) = S 1 x {x } x {0}. (5.20) 

Then V is a complete locally convex topological vector space. 
Take a (p e V expressed as 

<p(t,x,e) = (t,p(t,x,B),e) (5.21) 

satisfying 

D x /3(t,x Q ,0) = id, D n j3(t, x ,0) = (n>2), (5.22) 
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where j3 is a C w map. Define the representation n of K in V to be 

n(k)(y) = (ai(k) x D x a 2 (k)(t,x ,0) x id) o v o A(ky\ (5.23) 
and define Tp e V to be 

Tp = av(n)(<p) = f (a Y (k) x D x a 2 (k)(t, x Q , 0) x id) o ^ o A{ky l dk. (5.24) 

Then Thm.5.4 implies that Tp = n(k)(Tp) and this proves Eq.(5.15). Eqs.(5.16),(5.17) and 
(5.19) immediately follow from the definition (5.24) of Tp with Eqs.(5. 14),(5. 18),(5.20),(5.21) 
and (5.22). Since DTp(t, x ,0) = id, Tp is a C w diffeomorphism on a neighborhood of 
S 1 x {xq\ x {0} by virtue of the inverse mapping theorem. ■ 

Proof of Thm.5.1. Suppose that an infinite order RG transformation a t (y) is analytic 
in U. Then, so is an infinite order RG equation. Since the RG equation is invariant un- 
der the translation of t, Eq.(5.1) is invariant under the action defined by pulling back the 
translation by the RG transformation: 

{t, x) i-> (t + k, a t+ k o a~ l (x)) 

= (t + k,x + s(u^ k (x) - uf\x)) + 0(s 2 )), k£R. (5.25) 

Since a t is T-periodic, this defines the T 1 action on S 1 x M of the form of (5.5). 

Conversely, suppose that Eq.(5. 1) is invariant under the T 1 action (5.5). Let us rewrite 
Eq.(5.1) as 



V 




1 




X 




£g(s, X, S) 


(5.26) 






{ o J 





Then, the action (5.5) induces the action on S 1 x M x / of the form 

A{k) : (s, x, s) i-> (s + k, a 2 (k)(s, x, e), e), keR, (5.27) 

where a 2 (k)(s, x, e) = x + scr k (s, x, s). Since a 2 (k) satisfies Eqs.(5.14) and (5.18) for any 
x 6 M, Thm.5.5 applies to show that there exist an open neighborhood U = U Xo of 
S 1 x {x } x {0} and a C w diffeomorphism Tp satisfying Eqs.(5.15) to (5.17) and (5.19) for 
each xq 6 M. 

By taking a local coordinate near xq, we put x = and identify a neighborhood of 
xq with a neighborhood of in T XQ M. Then, Eqs.(5.16), (5.17) and (5.19) prove that Tp is 
expressed as 

Tp(t, x, s) = (t, x + sa.T,(t, x, s), e), (5.28) 
where a 3 is a C w map. Thus Tp defines the C w transformation \J/ t by 

x = \ff t (y) = y + sa 3 (t,y, e). (5.29) 
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Now Eq.(5.15) proves that if we transform Eq.(5.1) by x = if/ t (y,s), then the resultant 
equation is invariant under the action 

a\(k) x D x a 2 (k)(t, x Q , 0) x id : (t, x, s) i-» (t + k, x, s). (5.30) 

Since this is the translation of t, the resultant equation has to be an autonomous system 
of the form 

%=efR(y,s), (5.31) 
at 

where K is analytic in y and s. 

Finally, we show that there exist an RG transformation and an RG equation which co- 
incide with Eqs.(5.29) and (5.31), respectively. Let Eqs.(5.3) and (5.4) be one of the RG 
transformations and the RG equations for Eq.(5.1), respectively. Since both of Eqs.(5.31) 
and (5.4) are obtained from Eq.(5.1) by C°° transformations tff t and a t , respectively, there 
exists an C°° transformation cf> which brings Eq.(5.31) into Eq.(5.4): 

a, o <p(y, s) = iptiy, s). (5.32) 

Because of Claim 2.16, there exists an RG transformation c? f other than a, such that 
a, = a t o <p. This proves that there exists an analytic RG transformation a t = if/,. ■ 



5.3 Infinite order RG equation for linear systems 

In this subsection, we consider an n-dimensional linear system of the form 

x = sA(t, s)x, , 



eAi(t)x + s A 2 (t)x + 



x e C", 



(5.33) 



where A(t, s) is an nxn matrix. If this system satisfies the assumption (A), the RG method 
is applicable. Since the RG equation and the RG transformation are linear in y e C", it is 
convenient to define matrices Rj and up to be 



1 r' 

Ri = lim - Ai(s)ds, 

t^oo t J 

u { ? = j (Ai^-R^ds, 



(5.34) 
(5.35) 



and 



R t = lim 

t— >oo 



i rH 1 1 

- IJ^A^uf +A i (s)-J]uf ) R^ k 
J \k=i k=\ ) 

/ t r i-i ;-i 
J] Ai_ k (s)uf + Ms) - ^ ufRi- k - R 
<k=i k=l 



ds, 



ds, 



(5.36) 
(5.37) 
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for i = 2, 3, • • • , respectively. With these matrices, the m-th order RG equation and the 
m-th order RG transformation for Eq.(5.33) are defined by 



y = sRiy + s 2 R 2 y + • • • + s m R m y, yeC n , (5.38) 
a^y = y + eu^y + ■■■ + s m u ( t m) y, yeC", (5.39) 

respectively. Since the RG equation is a linear system with a constant coefficient, we can 
easily determine the stability of the trivial solution x = of Eq.(5.33) by using the RG 
equation. See Chiba [11] for the detail. 

In what follows, we suppose the following assumption : 

(L) The matrix A(t, s) is T-periodic in t and analytic with respect to s e D, where DcC 
is an open neighborhood of the origin. 

Let us consider the infinite order RG equation y = ^T=i £k Rky '■= f^(.e)y and the 
infinite order RG transformation a t y = y + J^ =i e^uf^y. Then, the fundamental matrix 
X(t, s) of Eq.(5.33) is given by 

X(t,s) = a t -e nE)t . (5.40) 

An initial value is given by X(0, s) = a = id + J^ =i e^u^ and it can be taken arbitrarily 
by choosing integral constants in Eqs.(5.35) and (5.37) appropriately. We suppose that 
the integral constants are chosen so that an initial value a is analytic in e e D. Then, 
X(t, s) is analytic in s e D. Now the question arises whether a t and K{e) are analytic in 
s. 

Theorem 5.6. Suppose that Eq.(5.33) satisfies the assumption (L) and ao is chosen to be 
analytic ins e D. Then, there exists a positive number r such that K{e), e^' and a t are 
analytic on the disk \e\ < r . 



Proof. Since a t is T-periodic (Lemma 2.1 (ii)), the equality 



X(t + T,e) = a t+T ■ e 

JR(£)t . JR(e)T 



holds. Putting t = yields 



= X{t,s)-e ns)T (5.41) 



e K(s)T = X(Qj e) -\ . X( j^ g) (5 42) 



Since X(t, s) is nonsingular and analytic in e e D for all t e R, is also analytic in 
s e D. Since det e n( - £)T + 0, the theory of analytic matrix-functions concludes that %{s)T 
is analytic on some disk {s e C | \e\ < r } c D (see Yakubovich and Starzhinskii [52], 
Erugin [17]). This proves that < R(e),e R ^ t and a t = X(t,s)e~' R( -^ t are also analytic on the 
disk. ■ 
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A few remarks are in order. Thm.5.6 with Thm.5.1 shows that Eq.(5.33) with the 
assumption (L) admits a Lie group action other than the scalar multiple : x h-» kx, k e C. 

In general, analyticity of e K( - £)T on D does not conclude analyticity of e n{£)t on D for 
all (eR. For example, consider the case K{e) = log(l + e), T = 1. 

The convergence radius r of !R(e) is given as follows: Fix so & D and a pass / in D 
from the origin to s . Let • • • , /t„(e) be eigenvalues of e^ (e)r . Suppose that there are 
z j such that /t ; (eo) = Aj(so) and passes of Aj(s) and /l/e) along / surround the origin 
(see Fig. 8). Then, log A,(eo) and log/l/eo) are located in different sheets of the Riemann 
surface. The smallest absolute value r = |eol among such s 's gives the convergence 
radius. 






logXj 


o 1 


log^ 



Fig. 8: Passes of eigenvalues /!,, Aj on C and passes of log A u log/l, on the Riemann 
surface. 



Floquet theorem states that for a given linear system with a periodic coefficient, there 
exist a periodic matrix Q{t) and a constant matrix B such that the fundamental matrix 
X(t) is written as X(t) = Q{t)e Bt . The RG method just gives the matrix B; B = Hie) = 
YX=\ £k Rk- Since Q(t) = a t is periodic, the stability of the trivial solution x = is 
determined by eigenvalues of B = TI{e), called Floquet exponents. 

Let T = 2n for simplicity. The matrix A(t, s) is expanded in a Fourier series as 
Mf,e) = 2-M C n( £ ) e '"'- By changing the independent variables as z = e mt , Eq.(5.33) is 
written as 

, +oo 

-^ = -wY J c n (s)z n ' l x. (5.43) 

— oo 

This is a linear system on a complex domain having the singularity at the origin. It is 
known that the fundamental matrix of this system is expressed as 



X(z) = S(z)-e 



M-logz 



(5.44) 



where S(z) is a single- valued matrix and M = M{e) is a constant matrix called mon- 
odromy matrix. It is easy to verify that ^(e) = iM and the RG method provides the 
monodromy matrix. 
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